Appendix A. 


Final Source Listing of Fuzzy Learning Modules 

for 

Shuttle Translational Controller 







SIMULATION APPLICATION: ARIC Translational Controller Simulation 


Phi_dot vs TIME 

RUN: V Bar Approach 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


2 


Tbc Nov 17 1992 02:53:45 PM 


SIMULATION APPLICATION: ARIC Ttauislational Controller Simulation 


push vs TIME 

RUN: V Bar Approach 



0 200 400 600 800 1000 1200 1400 1600 1800 

TIME (sec) 


MODULE: ORBITERJjn_rtng e 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


3 


Tic Nov 17 1992 02:53:45 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


failure vs TIME 

RUN: V Bar Approach 



TIME (sec) 

MODULE: ORBITER-lnwinge 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


4 


Ibe Nov 17 1992 02:53:45 PM 


SIMULATION APPLICATION: ARIC Thmslational Controller Simulation 


x[0] VS TIME 

RUN: V Bar Approach 



TIME (sec) 

MODULE: ORBl 1 hKlm_r*nge 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


5 


Tie Nov 17 1992 02:53:45 PM 



400 600 800 1000 1200 1400 

TIME (sec) 


MODULE: ORBlI , hKim_range 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


6 


Tbe Nov 17 1992 02:53:45 PM 






SIMULATION APPLICATION: ARIC Translational Controller Simulation 


v vs TIME 

RUN: V Bar Approach 



MODULE: ORBITER-lm_rangc 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


8 


Tbe Nov 17 1992 02:53:45 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[12] vs TIME 

RUN: V Bar Approach 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


Tie Nov 17 1992 02:53:45 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[25] vs TIME 

RUN: R Bar Approach 



TIME (see) 

MODULE: ORBITERini_mige 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


13 


Tbe Nov 17 1992 03:06:49 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[36] vs TIME 

RUN: R Bar Approach 



TIME (sec) 

MODULE: ORBITERJni_i»n*e 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


14 


Tie Nov 17 1992 03:06:49 PM 






SIMULATION APPLICATION: ARIC Thmslational Controller Simulation 


d[48] vs TIME 

RUN: R Bar Approach 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


16 


Tie Nov 17 1992 03:06:49 PM 





SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[49] vs TIME 

RUN: R Bar Approach 



TIME (sec) 

MODULE: ORBi 1 bKim_ringc 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


17 


Tic Nov 17 1992 03:06:49 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


unusualness vs TIME 

RUN: R Bar Approach 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


sfuel vs TIME 

RUN: R Bar Approach 



DATA SAMPLING FREQUENCY: 0200 Hi 


ORBnAL OPERATIONS SIMULATOR 


Tie Nov 17 1992 03:05:15 PM 


B3. Fly-Around from V-bar to R-bar vector in 30 minutes 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH X POS vs TIME 

RUN: Fly Around - V Bar To -R Bar 



VEHICLE: ORBITER.iute 
TARGET VEHICLE: SOLMAX.sUte 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


1 


Fri Nov 20 1992 02:23:36 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Y POS vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TARGET VEHICLE: SOLMAX.stite 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


2 


Fri Nov 20 1992 02:23:36 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Z POS vs TIME 

RUN: Fly Around - V Bar To -R Bar 


CURV 

LVLH 

ZPOS 

(feet) 



VEHICLE: ORBITER.sute 
TARGET VEHICLE: SOLMAX.jute 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


3 


Fri Nov 20 1992 02:23:36 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH X VEL vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TARGET VEHICLE: SOLMAX.iute 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


4 


Fri Nov 20 1992 02:23:36 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Y VEL vs TIME 

RUN: Fly Around - V Bar To -R Bar 


CURV 

LVLH o.OOl 
Y VEL 
(ft/s) 

0 



10.5 14 175 21 245 28 


TIME (minutes) 


VEHICLE: ORB ITER. «t»ie 
TARGET VEHICLE: SOLMAX.sute 
DATA SAMPLING FREQUENCY: 0.521 Hz 


Fri Nov 20 1992 02:23:36 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Z VEL vs TIME 

RUN: Fly Around - V Bar To -R Bar 


0.25 

0.2 

0.15 

0.1 

0.05 

-6.93889c- 17 

-0.05 

CURV _o l 

LVLH -015 
ZVEL 

(ft/s) -0-2 

-0.25 
-0.3 
-0.35 
-0.4 
-0.45 
-0.5 
-0.55 
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3 -3 7 10.5 14 17.5 21 24.5 28 

TIME (minutes) 


VEHICLE: ORBITER.sute 
TARGET VEHICLE: SOLMAX.sute 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


6 


Fri Nov 20 1992 02:23:36 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Y POS vs X POS 

RUN: Fly Around - V Bar To -R Bar 



VEHICLE: ORBITER.sute 
TARGET VEHICLE: SOLMAX.sute 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


7 


Fri Nov 20 1992 02:23:36 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Z POS vs X POS 

RUN: Fly Around - V Bar To -R Bar 



TARGET VEHICLE: SOLMAX.iute 
DATA SAMPLING FREQUENCY. 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


8 


Fri Nov 20 1992 02:23:36 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Y VEL vs X VEL 

RUN: Fly Around - V Bar To -R Bar 



VEHICLE: ORBITER.sute 
TARGET VEHICLE: SOLMAX.sute 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


9 


Fri Nov 20 1992 02:23:36 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Z VEL vs X VEL 

RUN: Fly Around - V Bar To -R Bar 



-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 -2.77556e-17 0.1 0.2 0.3 0 

CURV LVLH X VEL (ft/s) 


VEHICLE: ORB ITER. it»te 
TARGET VEHICLE: SOLMAX.«ute 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


10 


Fri Nov 20 1992 02:23:36 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


LVLH EULER PYR ROLL vs TIME 

RUN: Fly Around - V Bar To -R Bar 



0 3.5 7 10.5 14 17.5 21 24.5 28 31.5 

TIME (minutes) 


VEHICLE: ORB ITER, state 

DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


1 


Fri Nov 20 1992 02:23:55 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


LVLH EULER PYR PITCH vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.521 Hz 


Fri 


Nov 20 1992 02:23:55 PM 


ORBITAL OPERATIONS SIMULATOR 


2 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


LVLH EULER PYR YAW vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (minutes) 

VEHICLE: ORB ITER, state 

DATA SAMPLING FREQUENCY: 0.521 Hi 


ORBITAL OPERATIONS SIMULATOR 


3 


Fri Nov 20 1992 02:23:55 PM 




SIMULATION APPLICATION: ARIC Translational Controller Sim ulation 


LVLH BIASED BODY ROLL RATE vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:23:55 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


LVLH BIASED BODY PITCH RATE vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (minutes) 


VEHICLE: ORB ITER, stale 

DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


5 


Fri Nov 20 1992 02:23:55 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


LVLH BIASED BODY YAW RATE vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


6 


Fri Nov 20 1992 02:23:55 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


sfuel vs TIME 

RUN: Fly Around - V Bar To -R Bar 



MODULE: ORBl 1 hK. primary 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:24: 10 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


Phi vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


1 


Fn Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


Phi_dot vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


2 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


push vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:25:02 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


unusualness vs TIME 

RUN: Fly Around - V Bar To -R Bar 


0.0002 


0.00018 


0.00016 


0.00014 


0.00012 





200 400 


800 1000 1200 1400 1600 


TIME (sec) 


MODULE: ORBITER.lm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:25:02 PM 



SIMULA7TON APPLICATION: ARIC TiamJalionaJ Conm,ll CT S imulation 


failure vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITE3Um_clev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


5 


Fri Nov 20 1992 02:25:02 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


x[0] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


6 


Fri Nov 20 1992 02:25:02 PM 






8 











SIMULATION APPLICATION: ARIC Translational Controller Sim ulati on 


f[8] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:25:02 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[9] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



MODULE: ORBITERimelev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


14 


Fri Nov 20 1992 02:25:02 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[10] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


15 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[13] vs TIME 

RUN: Ry Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBllhK.lm_cIcv 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


16 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational 


Controller Simulation 


RTTM VS ™E 

RUN: Fly Around - V Bar To -R Bar 


0.990064 » 



MODULE: ORBITERJni_dev 

DATA SAMPLING Frequency.. 0.200 Hz 


00 1000 1200 — — 1 - 

14 00 i 6 oo 

TIME (see) 


ORBITAL OPERATIONS SIMULATOR 


17 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[ 15] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER-lm_clev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


18 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Thimlatlonal Commllar Simulauor 


d[9] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE. ORBITERJm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


19 


Fri Nov 20 1992 02:25:02 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[ 10] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITERIm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


20 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[12] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITERJra^elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


21 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[13] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORB ITERim _el cv 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


22 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[15] VS TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


23 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[16] vsTTME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_clcv 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


24 


Fri Nov 20 1992 02:25:02 PM 




0.981091 


MODULE: ORBITERJni_dev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


800 *°°° ™ i IwTST 

TIME (sec) 


ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:25:02 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[25] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 


MODULE: ORBlTEKJm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


26 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[27] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITERJm.clcv 

DATA SAMPLING FREQUENCY: 0.200 Hi 


ORBITAL OPERATIONS SIMULATOR 


27 


Fri Nov 20 1992 02:25:02 PM 





SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[28] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


28 


Fri Nov 20 1992 02:25:02 PM 







SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[31] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_clcv 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


30 


Fri Nov 20 1992 02:25:02 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[39] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


31 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[40] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


32 


Fri Nov 20 1992 02:25:02 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[42] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 


MODULE: ORBITERlm_clev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


33 


Fri Nov 20 1992 02:25:02 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[43] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (see) 

MODULE: ORB ITER. lm_elcv 

DATA SAMPLING FREQUENCY: 0.200 Hz 


OR3ITAL OPERATIONS SIMULATOR 


34 


Fri Nov 20 1992 02:25:02 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[45] VS TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


35 


Fri Nov 20 1992 02:25:02 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


0.979338 


d[46] vs TIME 

RUN: Fly Around - V Bar To -R Bar 


aflr&J383 



0.975428 


MODULE: ORBITER.lm_elcv 

DATA SAMPLING FREQUENCY: 0.200 Hz 


800 1000 1200 1400 1600 1800 

TIME (sec) 


ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:25:02 PM 








400 600 800 1000 1200 1400 

TIME (sec) 


MODULE: ORB ITER. lmrange 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


2 


Fn Nov 20 1 992 02:25: 1 1 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


push vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


3 


Fri Nov 20 1992 02:25: 11 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


failure vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORB ITER. lm_range 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


4 


Fri Nov 20 1992 02:25:1 1 PM 



SIMULATION APPLICATION: ARIC TrapsIaHonaJ Conpolle, S .muUuion 


x[0] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



ORBITAL OPERATIONS SIMULATOR 


5 


Fri Nov 20 1992 02:25: 1 1 PM 







SIMULATION APPLICATION: ARIC Translational Controller Simulation 


v vs TIME 

RUN: Fly Around - V Bar To -R Bar 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[8] VS TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (see) 

MODULE: ORBITERim^imge 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


9 


Fri Nov 20 1992 02:25: 1 1 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[12] vs TTME 

RUN: Fly Around - V Bar To -R Bar 



TTME (sec) 


MODULE: ORBITERJm_rangc 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


10 


Fri Nov 20 1992 02:25:11 PM 





SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[24] vs TIME 

RUN: By Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_rmnge 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


12 


Fri Nov 20 1992 02:25 : 1 1 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[25] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 


MODULE: ORBi l hK Jm_nmge 

DATA SAMPUNG FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


13 


Fri Nov 20 1992 02:25:11 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[36] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



MODULE: ORBITER.lm_range 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


14 


Fri Nov 20 1992 02:25: 1 1 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


0.962749 


d[37] vs TIME 

RUN: Fly Around - V Bar To -R Bar 


6PGJ827 



0.958905 


MODULE: ORBITERlm_range 

DATA SAMPLING FREQUENCY: 0.200 Hz 


800 1000 1200 1400 1600 1800 


TIME (sec) 


ORBITAL OPERATIONS SIMULATOR 


Fn Nov 20 1992 02:25:11 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[48] vs TIME 

RUN: Fly Around - V Bar To -R Bar 








TIME (sec) 

MODULE: ORBITERJm_range 

DATA SAMPLING FREQUENCY: 0.200 Hz 


Fri Nov 20 1992 02:25:1 1 PM 


ORBITAL OPERATIONS SIMULATOR 


17 




TIME (sec) 


MODULE: ORB ITER. lm_range 

DATA SAMPLING FREQUENCY: 0.200 Hz 



ORBITAL OPERATIONS SIMULATOR 


18 


Fri Nov 20 1992 02:25:1 1 PM 




TIME (sec) 


MODULE: ORBITER.lm_azim 

DATA SAMPLING FREQUENCY: 0.200 Hz 




ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:24:23 PM 





SIMULATION APPLICATION: ARIC Translational Controller Simulation 


unusualness vs TIME 

RUN: Fly Around - V Bar To -R Bar 



0 200 400 600 800 1000 1200 1400 1600 1800 

TIME (sec) 


MODULE: ORBITER.lm_tzim 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


4 


Fri Nov 20 1992 02:24:23 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


failure vs TIME 

RUN: Fly Around - V Bar To -R Bar 



DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


5 


Fri Nov 20 1992 02:24:23 PM 



0 200 400 600 800 1000 12 

TIME (sec) 


MODULE: ORBITER.lm_azim 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


6 









0 200 400 600 800 1000 1200 1400 1600 1800 

TIME (sec) 

MODULE: ORBITER.lm_azim 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 8 Fri Nov 20 1 992 02:24:23 PM 



0.015 



TIME (sec) 

MODULE: ORBITER.lm_azim 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


Fri Nov 20 1992 02:24:23 PM 










MODULE: ORBnrnim,^ TIME (sec) 

DATA sampling FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


11 


FnNov 20 1992 02:24:23 PM 




SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[8] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_azim 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


12 


Fri Nov 20 1992 02:24:23 PM 





SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[13] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lm_azirn 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


14 


Fri Nov 20 1992 02:24:23 PM 



TIME (sec) 

MODULE: ORBITER.lm_*zim 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


15 


Fh Nov 20 1992 02:24:23 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[ 10] vs TIME 

RUN: Fly Around - V Bar To -R Bar 



TIME (sec) 

MODULE: ORBITER.lni_azirn 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


16 


Fri Nov 20 1992 02:24:23 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC TVanslational Contralto Simulation 
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B4. Station-Keeping for 30 minutes at 50 feet distance on V-bar 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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CURV 
LVLH 
X POS 
(feet) 


31 


TIME (minutes) 



VEHICLE: ORBITER.sute 
TARGET VEHICLE: SOLMAX.tuie 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


Sat Nov 21 1992 04:22:37 PM 


SIMULATION APPLICATION: ASIC Tnmsiadonal Conner SimulaUon 


TARGET CENTERED ROTATING CURVILINEAR 

RUN: Station Keep At 200 Feet 


LVLH Y POS vs TIME 



TARGET VEHICLE: SOLMAX.statc 
DATA SAMPLING FREQUENCY: 0-521 Hz 


ORBITAL OPERATIONS SIMULATOR 


2 


Sat Nov 21 1992 04:22:37 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 


TARGET CENTERED ROTATING CURVILINEAR LVLH Z POS vs TIME 

RUN: Station Keep At 200 Feet 



VEHICLE: ORB ITER, stale 
TARGET VEHICLE: SOLMAX.sUte 
DATA SAMPLING FREQUENCY: 0.521 Hz 


ORBITAL OPERATIONS SIMULATOR 


3 


Sat Nov 21 1992 04:22:37 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


LVLH BIASED BODY ROLL RATE vs TIME 

RUN: Station Keep At 200 Feet 


2.03288e-19 

-5e-05 

- 0.0001 

-0.00015 

- 0.0002 

-0.00025 

LVLH 

BIASED-0-0003 













— 































k 















fSw 



BODY_o 00035 
ROLL 

RATE -0.0004 
(deg/s) -0.00045 
-0.0005 
-0.00055 




— 




— 










^ V 

































-w 














- 

-0.0006 ■ 














-0.00065 - 








J 



- 

~ 


-0.0007 - 











X 

— 

- 

n 

T < 

5 ^ 


« a 







14 17.5 21 

TIME (minutes) 


VEHICLE: ORBITER.sute 

DATA SAMPLING FREQUENCY: 0.52! Hz 


ORBITAL OPERATIONS SIMULATOR 


4 


S« Nov 21 1992 04:22:22 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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IMULATION APPLICATION: ARIC Translational Controller Simulation 
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Range ARIC Learning Parameters - Inputs 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[37] vs TIME 

RUN: V Bar Approach 



TIME (sec) 

MODULE: ORBiihKJm_range 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


15 


Tie Nov 17 1992 02:53:45 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


unusualness vs TIME 

RUN: V Bar Approach 
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B2. Shuttle R-bar Approach from 400 feet to 50 feet 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Thmslational Controller Simulation 
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SIMULATION APPLICATION: ARIC Uaodtitional Controller Sf-.*. 


f[13] vsTTME 

RUN: R Bar Approach 



MODULE: ORBITERJni_Minj 


TIME (sec) 


DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


14 


The Nov 17 1992 03:05:42 PM 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[24] vs TIME 

RUN: R Bar Approach 



TIME (sec) 

MODULE: ORBITERJni_tajn 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


Tie Nov 17 1992 03:05:42 PM 



SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Thmslational Controller Simulation 
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RUN: R Bar Approach 
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SIMULATION APPLICATION: ARIC Thmslational Controller Simulation 
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RUN: R Bar Approach 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC TYanslational Controller Simulation 
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RUN: R Bar Approach 
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SIMULATION APPLICATION: ARIC Ttanslational Controller Simulation 


f[13] VS TIME 

RUN: R Bar Approach 
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SIMULATION APPLICATION: ARIC Translational 
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OTON APPLICATION: ARIC Translational Controller Simulation 





SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[12] vs TIME 

RUN: R Bar Approach 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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RUN: R Bar Approach 
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SIMULATION APPLICATION: ARIC Thmslational Controller Simulation 
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HON APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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[: ARIC Translational Controller Simulation 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC TVanslational Controller Simulation 


Phi_dot vs TIME 

RUN: R Bar Approach 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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RUN: R Bar Approach 
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SIMULATION APPLICATION: ARIC Thmslational Controller Simulation 
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1.0 Introduction : 


As part of the RICIS project # AR.06 activity, the reinforcement learning techniques developed 
at Ames Research Center are being applied to proximity and docking operations using the 
Shuttle and Solar Maximum Mission (SMM) satellite simulation. In utilizing these fuzzy 
learning techniques, we also use the Approximate Reasoning based Intelligent Control (ARIC) 
architecture, and so we use two terms interchangeable to imply the same. This activity is 
carried out in the Software Technology Laboratory utilizing the Orbital Operations Simulator 
(OOS). 

This report is the deliverable D3 in our project activity and provides the test results of the fuzzy 
learning translational controller. This report is organized in six sections. Based on our 
experience and analysis with the attitude controller, we have modified the basic configuration 
of the reinforcement learning algorithm in ARIC as described in section 2. The shuttle 
translational controller and its implementation in fuzzy learning architecture is described in 
section 3. Two test cases that we have performed are described in section 4. Our results and 
conclusions are discussed in section 5, and section 6 provides future plans and summary for 
the project 


2.0 Fuzzy Learning System Configuration 

In this section, we have described the changes/modifications we have made in the fuzzy 
learning algorithms within the ARIC framework. These changes were determined necessary 
based on our study so far to properly utilize the fuzzy learning techniques for space operations. 
The ARIC algorithm with these changes provides us a baseline that we will use for all 
translational modules - range, azimuth and elevation control. 

a. Remove bias from all inputs and rules. 

We have learned from the attitude controller performance that the bias term really limits the 
fuzzy controller's ability to perform, especially when the bias term is kept constant at 0.5 value 
through out the learning cycle. The bias was zeroed out in several test cases and results showed 
good improvements, particularly in learning the environment. Therefore, we have decided to 
remove bias from our future implementations. We have also discussed this issue with Dr. 
Berenji at Ames and he also agrees that if we do not need bias, then, it would be better for the 
learning system. 

b. Updates of d's and f s both depend on the belief value of premise as well as consequent part 
of the rule. 

In our earlier experiments, we observed that the variations in all d’s and all f s were the same, 
meaning that if d(l) changes by an amount, then d(2) will change by the same amount, even 
though, the initial values of d(l) and d(2) are different. This behavior is not acceptable in the 
sense that the coefficients d’s and f s should have different variations based on which rules 
fire. We examined the formulas that update these weights, and concluded that the best way to 
incorporate effects of rule firing is to include the strength in the update of coefficient f. 
Similarly, the d's should be updated using the belief value from the fuzzification of the premise 
side of the rule. Both of these changes were discussed in detail with Dr. Berenji, Dr. Lea and 
our team. It is very appropriate to utilize these changes in the algorithm, and therefore, we 
implemented these modifications. 


c. Crisp failure with protection. 


Our current implementation uses what we call crisp failure, e.g. 0 or -1. at 1.4 DB, but we 
must protect the weight updating from too much punishing. Since failure is based on a 
parameter value, it can continue to be a failure for many cycles. If the weights are updated for 
few cycles with a large punishment, the controller goes unstable and can not recover at any 
time during the mission. This is not allowed and will not be allowed for space operations. So 
we need to protect against too much punishing that can result in total catastrophe. 

d. Action changed to no action rather than reverse action. 

In the earlier version of ARIC, the Stochastic Action Modifier (SAM) changed the 'push' from 
one way to the other way, if the random probability exceeded the measure of confidence p 
calculated using the Action Selection Network coefficients. This means that the controller 
requires a positive torque, but the random probability at that time is larger than the measure of 
confidence p (output of Action Selection network), then SAM changed the positive torque to 
negative torque. { 'push' is changed to ’-push’ ). This process for space operations is not 
acceptable because the negative torque in such case will perturb the state too much and the 
controller may not be able to recover from such an action. 

In space operations, typically, an action is taken or no action is taken. It is seldom the case that 
an action is changed to its reverse. We also understand that such a process is desired in 
learning process, because the fuzzy learning technique is exploring all possible solution space. 
However, if we want to implement these learning techniques on a spacecraft so that it learns 
real time, this process must be changed to lower the risks. Therefore, we suggest to use the 
following procedure. 

If the random probability is larger than the measure of confidence p, then SAM should change 
the action to no action, meaning that the push value is set to zero, and not to its negative. We 
accomplish several benefits from this type of process. 

First, the controller will still learn from the search of solution. Now it is learning if the action 
recommended by the controller is really working or not It really leams what happens when the 
action is not taken. Second, the controller will not be punished too much and thus it will be 
possible to maintain the performance rather than degrade it. Third, the solution will be slowly 
discarded rather than catastrophic failure. The weight updates during this time ma y result in a 
better solution.. Fourth, process like this one is acceptable for space operations without 
increasing the risk of failure. Operationally, it is never acceptable to change an action to its 
reverse until there is a very high confidence that the reverse action is really needed. In our 
technique, since we have low confidence in an action, it should not be interpreted as high 
confidence in reverse action. Thus, logically it also makes sense for caring out operations. 

e. No changes to computations for measure of confidence. 

The measure of confidence 'p* is computed as originally proposed rather than normalizing 
using the number of rules. During the analysis of attitude controller performance we suggested 
that the measure of confidence p should be normalized per rule so that it will not continue to 
saturate during the learning process. We performed several attitude controller tests with this 
normalization, and had observed good learning rate for the neural network. However, now our 
main concern stems from the fact that the decision to fire jets was not arrived at by firing just 
one rule. Entire rulebase was utilized in the process, so for that reason we should not normalize 
the p-value. One can also argue that all rules were not used in deriving the controller action. We 
should properly scale the p-value using the number of rules fired. If we know that only four 
rules are responsible for the decision, then, we can normalize using these four rules. However, 


it is not possible to find out which four rules are responsible. Thus, it is best not to normalize 
the measure of confidence. Furthermore, the d's and fs are being updated using the firing 
strength from the left hand side and right hand side of the rules anyway. Thus, the calculation 
of ’p’ will include effect of long term behavior. 

f. Input parameters normalized between 0.0 and 1 .0. 

We were first normalizing the input parameters phi error and phi_dot error between -1.0 and 
+1.0. The reinforcement learning algorithms at this time require that the input parameters be 
normalized between 0.0 and +1.0. We have changed the attitude controller to reflect this 
change. In our opinion, there are no guidelines in neural network algorithms how to normalize 
input and output parameters. It is a research area to generate guidelines which parameters 
should scaled between 0.0 and 1.0, and which one should be scaled between -1.0 and +1.0. 
For example, range has no negative values, and there is what so ever no interpretation when a 
sensor measure negative range except that the sensor is failed. Similarly, mass of an object is 
never negative. 

g. Overlapping Membership Functions. 

During the analysis of attitude control performance, we discovered that the fuzzy learning 
technique uses a lot more fuel compared to fuzzy controller alone. Our results showed that the 
fuzzy learning or ARIC used three to eight times the fuel used by fuzzy only control. Further 
analysis showed that the defuzzification by Tsukomoto's method does not allow to use 
triangular membership functions for output. As a result we had to use only one side of the 
triangle. When we used only one side of the triangle, it resulted in a large hysterisis in the jet 
firing command sequence. Once the jets are turned on to provide torque in a desired direction, 
they are not turned off as soon as the rate error is corrected. The jets remain on for additional 
three to four cycles providing increased rate in the reverse direction. The angle error then hits 
the other side of the deadband. The control system then fires jets to reduce the angle error that 
was really caused by the previous unnecessary jet firing. Net effect is larger fuel usage and 
more angular activity. 

We further analyzed the rate error membership functions, and decided to overlap them (as 
shown in Fig. 1) to reduce the net hysterisis. We performed several tests to find a proper 
overlap, and now we have reduced the hysterisis to no hysterisis. We have decided to use this 
attitude controller in all our future tests. The rulebase for the attitude controller is shown in 
Table I for completeness. The translational controller uses this new attitude controller with 
overlap in rate membership functions. However, we will not utilize such overlap in rate errors 
for translational controller because we need to analyze the performance without such overlap. 
At this time we do not know if a hysterisis exists in the translational control. 



Attitude Error (degrees) 



Attitude Rate Error (degrees per second) 



desired firing pulse level 


KEY: 

NB - Negative Big, NM - Negative Medium, NS - Negative Small, ZO - Zero, 
PS - Positive Small, PM - Positive Medium, PB - Positive Big 


Fig. 1 Attitude Controller Membership functions - 
Overlapping Rate Error MBFs 


Table I. Fuzzy Rulebase for Attitude control 



Angle Error 

NB NM NS ZO PS PM PB 

NB 

PM PM PS PM 

NM 

PM PM PS PM 

Rate NS 

PS PS PS PS 

Error ZO 

PS PS PS ZO NS NS NS 

PS 

NS NS NS NS 

PM 

NM NS NM NM 

PB 

NM NS NM NM 


KEY: 

NB - Negative Big, NM - Negative Medium, NS - Negative Small, ZO - Zero, 
PS - Positive Small, PM - Positive Medium, PB - Positive Big 
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Fig. 2 Definition of Relative Trajectory Control 



Table III. Rulebase for azimuth control 


AZIMUTH RATE ERROR 
NM NS ZO PS PM 



PM 



NM 

NM 

NM 

AZIMUTH 

PS 



NS 

NS 

NM 

ANGLE 

ZO 

PM 

PS 

ZO 

NS 

NM 

ERROR 

NS 

PM 

PS 

PS 




NM 

PM 

PM 

PM 




Table II. Rulebase for elevation control 


ELEVATION RATE ERROR 
NM NS ZO PS PM 

PM 

NM NM NM 

ELEV. PS 

ZO NS NS NM 

ANGLE ZO 

PM PS ZO NS NS 

ERROR NS 

PM PS PS ZO 

NM 

PM PM PM 


KEY: 

NB - Negative Big, NM - Negative Medium, NS - Negative Small, ZO - Zero, 
PS - Positive Small, PM - Positive Medium, PB - Positive Big 
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Fig. 5 Membership Functions for Range Control 



Table IV. Rulebase for range control 


RANGE RATE ERROR 

NB NM NS ZO PS PM PB 


PB 

PM 

PS 

RANGE ZO 
ERROR NS 
NM 
NB 


NB NB NB 
NM NB NB 
NS NM NB 
PM PS PS ZO NS NS NM 
PB PM PS 
PB PB PM 
PB PB PB 


KEY: 


NB - Negative Big NM - Negative Medium, NS - Negative Small, ZO - Zero, 
PS - Positive Small, PM - Positive Medium, PB - Positive Big 




4.0 Description of Test Cases 


For our translational controller we have performed four testcases. These testcases were also 
performed for the fuzzy logic based 6 DOF controller and we have the man-in-the-loop 
simulation data that provides us insight as to how the Shuttle crew is trained for performing 
these proximity operations. 

The first test case is the v-bar approach where the shuttle is initialized at 400 feet distance in the 
front of the satellite. The pitch attitude of the Shuttle is 90 degrees (as shown in Fig. 6) so that 
a crew member can see the satellite from the window or the Crew Optical Alignment System 
(COAS) for measuring the angles of the satellite relative to the line of sight In our case we are 
using a radar like sensor placed at its proper location for measurements. The Shuttle has a 
closing rate of 0.4 feet per second. Its mission time-line calls for maintaining certain closing 
rate as a function of distance. As it approaches 50 feet distance, the mode changes to Station- 
keeping during which this distance is maintained. It approximately takes 1400 seconds for the 
Shuttle to complete this transition from 400 feet to 50 feet. The second test case known as r-bar 
approach is similar to the first one. Instead of v-bar it approaches along r-bar. The Shuttle is 
initialized at 400 feet on r-bar which is the z-axis of the LVLH frame as shown in Fig. 6. The 
pitch attitude of the Shuttle is zero degrees, and it sees the satellite from the window or COAS 
or radar. The trajectory is as shown in Fig. 6. 

The third test case is the Fly-Around test case. The Shuttle is initialized at 200 feet on v-bar 
with pitch attitude of 90 degrees very similar to the v-bar case. However, its pitch rate is 
initialized at 0.05 deg per sec, and thus it will rotate to increase the pitch angle. As the Shuttle 
rotates, the elevation angle will change and the range rate may also change. The elevation 
control rulebase will see this change in the elevation angle and fire the jets to move the Shuttle 
upward. As a results, the Shuttle will translate upward of v-bar. As the pitch rate continues to 
be non- zero, the elevation angle will continue to change, and the elevation control will keep the 
Shuttle pushing upward. When the Shuttle moves sufficiently upward, its range will increase, 
and then, the range control will initiate the range rate to reduce the range to 200 feet. Thus, the 
Shuttle will traverse a circular path as shown in Fig. 7, and will reach the negative r-bar in a 
certain given time. Actually, the initial pitch rate of 0.05 will set the arrival time of 1800 
seconds. If the 90 Fly-around is desired in 900 seconds, then, the initial pitch rate should be 
0.1 deg per sec. The mission time-line is typically set using this type of calculations. 

The fourth testcase is the station-keeping testcase, where the Shuttle is initialized at 200 feet 
distance on v-bar with 90 deg pitch angle and zero pitch rate. The range must be maintained 
within the range deadband as the mission time increases. Also, the elevation and azimuth 
angles must be maintained near zero. This activity is very simple and very boring for the 
Shuttle crew. There is not much jet firing activity and the drift in range is very slow, so the 
Shuttle is not moving relatively fast. However, this station-keeping is an important activity 
because it allows the Shuttle crew to check out sensors and other systems. We have set up this 
activity for 1800 seconds to be consistent with all other testcases. 





+r-bar 


Fig. 7 Definition of Proximity Operations fly-around Segment 


5.0 Results and Conclusions 


In all four testcases the translational controller in ARIC architecture works ; it achieves the 
desired position at the end of the proximity operations segment. This behavior is expected 
because the translational controller derived in ARIC is from the fuzzy logic based conffoller 
which has given very good results. In all testcases, we observed that the fuel usage is 
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this behavior is again expected for this controller. improvement. 
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6.0 Future Plans and Summary 
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Al. Source code for range control 







leam_range.c 

♦include <stdio.h> 

♦include <math.h> 

♦include <sys/types . h> 

/* EXTERNAL DATA STRUCTURE DEFINITION */ 
♦include . /orb_f uzzy/learn_cycle .h" 


♦define max(x,y) { (x 

♦define min(x,y) ( (x 

♦define Gamma 0 . 9 

♦define Beta 0.2 

♦define Beta_h 0.05 

♦define Rho 1-0 

/* 1 July 92 Change Rho__h 
♦define Rho_h 0.8 

♦define Rhol 2.0 

♦define Rho_hl 0 . 4 

extern double sgn(); 
extern double expo ; 
extern double rnd(); 
extern double sim_time{); 


y) 

y) ? 


? x 


: y ) 
y ) 


from 0.2 to 0.8 * / 


learn_range (L) 

LEAFN_CYCLE * L; /*IN : */ 

{ 

int i,il, j, k; 

double range_match ( ) , range_calculate_z_array () ; 


double temp ; 


L— >x [ 0 ] - <L->Phi+10) /20 ; 

L->x[l] - (L->Phi_dot + 5)/10.0 ; 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

/* 1 July 1992 - Set Bias to 0.0 */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

L->x[2] * 0.0 ; 

L->f ailure * 0; 

/* Set Up and evaluate the failure criteria */ 

if ( (f abs (L->Phi) >10) II (fabs (L->Phi_dot) > 5) ) { 

L->failure - -1.; 

if ( L->learn_flag ) { . . . 

fprintf (stderr, " learn_range : Failure at %f . \n", sim_time ( ) ) 

} /* end if */ 

} /* end if */ 

/ ★ CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 

/* CC 6 November, 1992 

/★ CC Turn off learning when there is a failure CC / 

/★ CC until there has been no failure for 4 pass CC */ 

/ * CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC * / 

L->fourth_pass_failure - L->third_j>ass_f ailure ; 

L->third _j>ass_f ailure « L->secondjpass_f ailure ; 

L->second_j5ass_f ailure “ L->f irst_pass_failure ; 

L->f irst_jpass_f ailure ® 0 ; 
if (L->failure " -1) { 

L->first_pass_f ailure - 1 ; 

} /* end if */ 

/* If No Failure in the last 4 passes then turn learning on */ 
if( !L->fourth_pass_f ailure && 


! L->third_pass_f ailure && 

! L->second_pass_f ailure £& 
! L->f irst_ - pass_f ailure ) { 

L->learn_f lag - 1 ; 

} /* end if */ 


LJangep 


/* If Failure first pass, but not second pass, let the learn 
flag stay on, so that the failure is processed. On the next 
pass (where there is a failure second pass, but not third 
pass) then turn learning off until there are no failures for 
four passes */ 

if { L->second__pass_f ailure £& ! (L->third_j>ass_f ailure) ) { 

L->learn_f lag - 0 ; 

) /* end if */ 


/* output: state evaluation */ 

for (i - 0; i < 25; i++) 

{ 

L->sum * 0.0; 

for ( j - 0; j < 3; j++) 

{ 

L->sum +- L->a [i*3+ j ] * L->x_old[j]; 

J 

L->y[i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; 

> 

L->sum - 0.0; 

for(i -0; i < 25; i++) 

{ 

L->sum +— L->c[i) * L->y[i]; 

) 

L->suml * 0.0; 

for ( j - 0;j < 3; j++) 

{ 

L->suml +« L->b[j] * L->x_old[j]; 

} 

L->v ” L->sum + L->suml; 

/* output: action */ 
for(i * 0; i < 25; i++) 

( 

il-i; 

L->w[i] — rangejmatch (il, L) ; 

L->zl[i] - range_ calculate_z_array (il, L) ; 

) 

L->numl - 0.0; 

L->denom - 0.0; 
ford “0; i < 25; i++) 

{ 

L->numl +- L->w [i] * L->zl[i] * L->f[i] ; 
L->denom +— L->w [i] *L->f [i] ; 


) 

/* JUNE 9, 1992 - CORRECTION ! ! ! */ 

/* Add test for denom very small compared to suml - no rule firing zone */ 
/* L->push - (1000 . 0*fabs (L->denom) <fabs (L->suml) ) ?0 . 0 :L->suml/L->denom; */ 
if (fabs (L->numl) <-0 . 01 I I fabs (L->denom) <-0 . 01) ( 

L->push -0.0 ; 

) else { 

L->push - L-> numl/L->denom ; 

) /* end if */ 



'^ieanucange.c 

/* output: action computations completed */ 
for (i - 0; i < 25; i++) 

{ 




L->sum - 0.0; 

for (j - 0; j < 3; j++) 

L->sum +- L->d[i*3+ j] * L->x_old[j]; 

L->z[i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; 

} 

L->sum2 - 0.0; 

L->sum3 ” 0.0; 

ford -0; i < 3; i++) 

L->sum2 +- L->e [i] * L->x_old[i]; 

for (i — 0;i < 25; i++) 

L->sum3 +- L->f [i] * L->z [i] ; 

/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


/* L->unusua lne s s = (L->push >0) ? 1.0 — L->p . P* 

/ * cccccccccccccccccccccccccccccccccccccccccccccccccccc * / 

/* 28 Sept 1992 - Set stocastic action to zero out cmd */ 

/ * cccccccccccccccccccccccccccccccccccccccccccccccccccc * / 

/* L->push - ( rnd() <» ( (L->p+l . 0) /2 . 0) ) ? L->push : -L->push; / 

/ * ccccccccccccccccccccccccccccccccccccccccccccccccccccc 

28 Sept 1992 - Change definition of unusualness to 
correspond with Hamid' s notes . 

L->push - ( rnd() <- ( (L->p+l . 0) /2 . 0) ) ? L->push : 0.0 ; 

L->unusualness - (L->push >0) ? 1.0 - L->p : -L->p, 

*/ 

if( md() <- (<L->p+1.0)/2.0) ) { 

L->unusualness *1.0 - L->p ; 

} else { 

L->push • 0.0 ; 

L->unusualness * -L->p ; 

} /* end if */ 


CC JULY l f 1992 - Change Normalize of sum4 
CC L->sum4 - L->sum3 + L->sum2; 

CC JULY l r 1992 - Normalize of sum3 by # rules 
CC L->sum4 * L->sum3 / 31.0 + L->sum2 / 3.0 ; 


/’ 


CC 

CC 

CC 

CC 


:c OCT 2, 1992 - Do NOT normalize for # rules - • - ^ 

ZCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/ 

»sum4 * L->sum3 + L->sum2 / 3.0 ; 

;,->p - l.o / (1.0 + exp (-1.0 * L->sum4) ) ; 

/ * CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC * / 

15 April 1992 - Use temp variable - not push */ 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 

t — snnoVi m i mH ( \ L— >d) ? L*->push i — L-’^’push; */ 


/* using new input values and unmodified weights. */ 

/* Use y_new and v_new so not to destroy y and v. */ 

for (i - 0; i < 31; i++> 

{ 

L->sum * 0.0; 

for ( j - 0; j < 3; j++) 




{ 

L->sum +- 

} 

L->y_new [ ij - 

} 


^angeic 


L->a [i*3+ j] * L->x [ j] ; 

1.0 / (1.0 + exp(-1.0 * L->sum) 


L->sum — 0.0; 
L->suml - 0.0; 
L->sum2 - 0.0; 



for ( j - 0; j < 3; j++) 

L->suml +■» L->b[j] * L->x [ j ] ; 

for(i - 0; i < 25; i++) 

L->sum2 +« L->c[i] * L->y_new[i]; 

L->sum - L->suml + L->sum2; 

L->v__new ■» L->sum; 


/ *ccccccccccccccccccccccccccccccccccccccccccccccccccccccccc* / 

/* 17 September, 1992 */ 

/* This logic depends upon a "crisp” (two valued) failure */ 

/* It has been replaced by "fuzzy" failures */ 

/* action evaluation */ 
if (L->failure) 

L->r_hat - L->f ailure - L->v; 
else 

L->r_hat - L->failure + Gamma * L~>v_new - L->v; 

/ *ccccccccccccccccccccccccccccccccccccccccccccccccccccccccc * / 

/* L->r_hat - -1.0 * L->failure + 

f* ( 1.0 - L->f ailure ) * Gamma * L->v_new - L->v ; 

/* modification and update to parameters */ 

for(i - 0; i < 25; i++) 

( 

L->f actor 1 - Beta_h * L->r_hat * L->y[i] * (1.0 - L->y[i]) * sgn(L->c[ 
L— >c[i] +“ Beta * L->r_hat * L— >y[i]; 

for ( j =0; j < 3; j++) 

{ 

L->a[i*3+j] +- L->factorl * L->x_old[j]; 

} 

) 


for(i - 0; i < 3; i++) 

L->b[i] +- Beta * L->r_hat * L->x_old[i]; 

for(i - 0; i < 25; i++) 

{ 

/♦cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

CC 7 July, 1992 : Weight updates of D's by rule firing strength CC 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc*/ 

L->factor2 - Rho_h * L->r_hat * L->z[i] * (1.0 - L->z[i]) * sgn(L->f[i 

ualness ; 

for ( j - 0; j < 3; j++) { 
if < L->learn_f lag ) { 

!< - >d[i*3+j] +“ L->factor2 * L->x_old[j] * L->d_weights [i*3+ j] ; 

) /* end if */ 

) /* end for */ 

) 



for(i 


0; i < 25; i++) 
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< 



if( L->learn_f lag ) ( ^ T ^ 

L_>f[i] +— Rho * L->r__hat * L->unusualness L->z[i] 

} /* end if */ 


* L->w [ i] 


} 


for (i * 0; i < 3; i++) 


{ 


L->e[i] +* Rho * L->r_hat * L->unusualness * L->x_old[i] ; 


L->x_old[0] - L->x [ 0 ] ; 
L->x_old[l] - L->x[l] ; 
L->x old [2] - L->x [2] ; 








double sgn(x) 
double x; 

{ 

if (x < 0.0) 

return (-1.0); 
else if (x > 0 . 0) 
return (1.0) ; 
else 

return (0.0) ; 


/* zero one function returns 0 for negative numbers 

1 for values > 1 
x for values between 0 and 1 */ 

double zero_one(x) 


double x; 

if (x < 0) return (0.0); 
else if (x > 1) return (1.0); 
else return (x) ; 


I* ********** Membership Functions for Range (Phi) 

double range_nbl (x) 
double x; 

return (min ( max((-4-x)/4 , 0.0 ), 1.0 )); 

) 




double range_nml (x) 
double x; 

££ ( x <« -4.0) return (min ( max((x+8)/4 , 0.0 ), 

else return (min( max(( -x-2)/2 , 0.0 ), 1.0 )); 


) 

double range_nsl(x) 


1.0 )); 





if (x <- -2.0) return (min ( max ( (x+4.0)/2 , 0.0 ), 1 0 ))• 

else return (min< max ( -x/2 , 0.0 ), 1.0 )); 


double range_zol (x) 
double x; 

{ 

if (x <- 0) return (min ( max( (x+2)/2 , 0.0 ), 1.0 ))• 

else return (min ( max( (-x+2)/2 , 0.0 ), 1.0 )); 

double range_psl(x) 
double x; 

( 

if (x <- 2) return (min ( max ( x/2 , 0.0 ), 1.0 )); 

else return (min ( max( (-x+4)/2 , 0.0 ), 1.0 )); 

double range_pml(x) 
double x; 

( 

if (x <- 4) return (min ( max<( x-2)/2 , 0.0 ), 1.0 ))• 

else return (min ( max(( -x+8)/4 , 0.0), 1.0)); 


double range_pbl (x) 
double x; 

{ 


return (min ( max ( ( 

} 

/★*★★*★★*★********* 
double range__nb2 (x) 
double x; 

{ 


x-4) / 4 , 0.0 ), 1.0 )>; 

* Range-Rate Membership Functions 


return (min ( max ( (-0 . 5-x) / . 5 , 0.0 ), 1.0 )); 




*/ 


double range_nm2 (x) 
double x; 

{ 

if (x <- -.5) return (min ( max ( ( x+1.0)/.5 , 0.0 ), 1.0 )); 

else return (min ( max(( -x-.15)/.35 , 0.0 ), 1.0 )); 

double range_ns2 (x) 
double x; 

{ 

if (x <- -.15) return (min( max(( x+.5)/.35 , 0.0 ), 10))- 

else return (min ( max(( -x/.15) , 0.0 ), 1.0 )); 

double range_zo2 (x) 
double x; 

{ 

if (x <- 0) return (min ( max(( x+.15)/.15 , 0.0 ), 1.0 )); 

else return (min ( max(( -x+.15)/.15 , 0.0 ), 1.0 )); 

double range__ps2 (x) 
double x; 

/* Modified from x/.l to put gap between ns2 and ps2 */ 

/* 29 Sept 1992 */ 

if (x <- .15) return (min ( max (x/ . 15 , 0.0 ), 1.0 )); 

else return (min ( max(< -x+.5)/.35 , 0.0 ), 1.0 )); 

double range_pm2 (x) 
double x ; 
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if (x <- .5) return <min( max(( x-.15)/.35 , 0.0 ), 1.0 )); 

else return (min ( max ( ( - x+1 . 0 ) / . 5 , 0.0 ) , 1.0 )); 

} 




double range_j3b2 (x) 
double x; 

{ 

return (min ( max ( ( x-.5)/.5 , 0.0 ), 1.0 )); 

} 


/★***★★***★★★*★*★* Defuzzification Process with Delta-V Membership Functions 




double range_nb3 (x) 
double x; 

{ 

return (-0.2 - 0.3*x); 

} 


double range_nm3(x) 
double x; 

{ 

return (“0.1 - 0.1*x); 

} 

double range_ns3 (x) 
double x; 

( 

return (-0 . l*x) ; 

1 

double range_zo3 (x) 
double x; 

{ 

return (0.0) ; 

) 

double range_jps3 (x) 
double x; 

( 

return (0 . l*x) ; 

} 

double range_pm3 (x) 
double x; 

{ 

return (0.1 + 0.1*x); 

} 

double range_j>b3 (x) 
double x; 

( 

retum(0.2 + 0.3*x); 

) 


double range_match (i f L) 
int i; 

LEARN_CYCLE *L; /*IN : */ 

{ 

double temp; 
switch (i) { 


case 0: 




L->d_weights [i*3+0] - zero_one (range_nbl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero__one (range_zo2 (L->Phi_dot) *L->d [i*3+l] ) ; 

temp -min ( zero_one (range_nbl (L->Phi) *L->d[i*3+0] ) , zero_one (range zo2 (L->Phi dot) *3 
->d[i*3+l] ) ) ; 

return (temp) ; 

case 1 : ^ 

L->d_weights [i*3+0] - zero_ one (range_nbl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (range_ns2 (L->Phi_dot) *L->d [i*3+l] ) ; 
temp-mi n ( zero_one (range_nbl (L->Phi) *L->d [i*3+0] ) , zero_one (range ns2 (L->Phi dot)* 
>d [i*3+l] ) ) ; “ ^ 

return (temp) ; 
case 2 : 

L->d_weights [i*3+0] - zero^one (range_nbl (L->Phi) *L->d [i*3+0] ) ; 

L~>d_weights [i*3+l] - zero_one (range_nm2 (L->Phi_dot ) *L->d[i*3+l] ) ; ^ 

temp-min ( zero_one (range__nbl (L->Phi) *L->d[i*3+0] ) , zero_one (range_nm2 (L->Phi_dot) *T.- 
>d[i*3+l] ) ) ; 

return (temp) ; ^ 

case 3: 

L->d_weights [i*3+0] - zero_one (range_nbl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (range_nb2 (L->Phi_dot) *L->d [i*3+l] ) ; 
temp-min ( zero_one (range_nbl (L->Phi) *L->d [i*3+0] ) , zero^one (range nb2 (L->Phi dot)*** 
>d [i*3+l] ) ) ; 

return (temp) ; 
case 4: 

3»d_weights [i*3+0] - zero_one (range_nml (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (range_zo2 (L->Phi_dot ) *L->d [i*3+l] ) ; 
temp-min { zero^one (range_nml (L->Phi) *L->d[ 1*3+0] ) , zero_one (range zo2 (L->Phi dot)* 
>d[i*3+l] ) ) ; “ _ 

return (temp) / 
case 5 : 

L->d_weights [i*3+0] - zero__one (range_nml (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (range_ns2 (L->Phi_dot) *L->d[i*3+l] ) ; — * 

temp-min ( zero^one (range__nml (L->Phi) *L->d[i*3+0] ) , zero_one (range_ns2 (L->Phi dot) *L- 
>d [i*3+l] ) ) ; 

return (temp) ; 

case 6 : — 

L->d_weights [i*3+0] - zero_one (range_nml (L->Phi) *L->d[i*3+0) ) ; 

L->d — weights [i*3+l] - zero_one (range_nm2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero_one(range_nml (L->Phi) *L->d[i*3+0] ) , zero^one (range_nm2 (L->Phi_dot) *. 
>d [i*3+l] ) ) ; *“ 

return (temp) ; 
case 7 : 

L->d_weights [i*3+0] - zero_one (range_nml (L->Phi) *3»d[i*3+0] ) ; _ 

L->d_weights [i*3+l] - zero^one (range_nb2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero_one (range_nml (L->Phi) *L->d[i*3+0] ) , zero_one (range_nb2 (L->Phi_dot) *V 
>d[i*3+l] ) ) ; 

return (temp) ; w 

case 8 : 

L->d_weights [i*3+0] - zero_ one (range_nsl (L->Phi) *L->d(i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (range_zo2 (I»Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero_one (range_nsl (L->Phi) *L->d[i*3+0] ) , zero_o ne (range_zo2 (L->Phi_dot) * isr 
>d[i*3+l] ) ) ; 

return (temp) ; 
case 9 : 

L->d_weights [i*3+0] - zero_one (range_nsl (L->Phi) *L->d[ 1*3+0 ] ) ; 

L->d_weights [i*3+l] - zero_one (range_ns2 (L->Phi_dot) *L->d[i*3+l] ) ; 
temp-min ( zero_one (range_nsl (L->Phi) *L->d[i*3+0] ) , zero_one (range_ns2 (L->Phi_dot) *.: 
>d[i*3+l] ) ) ; _ 

return (temp) ; 
case 10: 

L->d_weights [i*3+0] - zero_one (range_nsl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (range_nm2 (L->Phi__dot) *L->d[i*3+l] ) ; — 

temp-min ( zero_one (range — nsl (L->Phi) *L->d[i*3+0] ) / zero_one (range__nm2 (L->Phi__dot) *L- 
>d[i*3+l] ) ) ; 
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return (temp) ; 

CaS L->d weights [i*3+0] - zero one (range_nsl (L->Phi) *L->dU*3+0] ) ; 

L_>d — weights [i* 3+1] - zero'one (range_nb2 <L->Phx_dot) *L->dU*3+l] ) . 

te^Snf zeio.one (range.nll (L->Phi) *L->d[i*3+0] ) , zero.one (range_nb2 <L- 

>d[i*3+l] ) ) ? 

return (temp) ; 

>d[i*3+l] ) ) ; 

return (temp) ; 

CaS L->diweights[i*3+0] - *ero_one (range jaljWPhij^-^U^+Ol^J ? 

T->d weicrhts [i*3 + l] *“ zero one { range_j?b2 (L->Phi_ dot ) ,, 

te^nf zero.one (range_pll (L->Phi) *L->d[i*3+0] ) , zero_one <range_j?b2 (L 

>d[i*3+l] ) ) ; 

return (temp) ; 

° aS L->d weights [i*3+0] - zero_one (range_ps! (L->Phi) *L-> d{i* ^+0]) ; 

L->d~ weights [i* 3+1] - zero one (range __pm2 (L->Phi_dot) L->dU ' _ .. 

teS=Snf zeio_one (rangej^l <L->Phi> -L->dli*3 + 0] ) , zero_one (range J»2 (L 

>d[i*3+l] ) ) ; 

return (temp) ; 

° aS L->d* weights [i*3+0] - zero one(range__psl (L->Phi) *L->dti*3+0] ) ; 

L->d — weights [i*3+l] - zero - one (rangej>s2 (L->Phi_dot) “L->d[i*3+l] , 

te2=£n! zero_one (r.nge_p7l (L->?hi) *l->d[i-3 + 0] ) . zero.one (range J>s2 (t 

>d[i*3+l] ) ) ; 

return (temp) ; 

CaS V>V weights [i*3+0] - zero.one <range_psl (L->Phi) J ” + i , , . 

temp-rrdn^zero^one^range^p^l^L-^Phit^L-^d^i^S+oT ) ? zero.one (range_zo 2 (L 

>d[i*3+l] ) ) ; 

return (temp) ; 

CaS L->d weights [i*3+0] - zero one (range_pml(L->Phi) *L->d[ r*3+0] ) , 

L L ->dIwei g ghts U‘3+1] ” ^ r p~°^->phit^->d^i*3+0T^? t zero > one (^>2 (h 
temp-min ( zero_one (range_pml (L->Pnij u "J" - 

>d[i*3+l] ) ) ; 

return (temp) ; 

CaS L->d^.weights [i*3+0] - zero.one < r.nge_pml (L->Phi) .L->dIi*3t0n + ; ^ ^ . 

. zerc.one .ranged « 

>d[i*3+l] ) ) ; 

return (ten^) ; 

CaS L->diweights [ i* 3+0 ] - zero_one j r »nge^ml jL->Phi) , , 

tl^nf zeio’oneulng.l^ia-'phi) 'I.->d!i*3 + 0] ) , s«ro_on. <range_ps2 

>d[i*3+l] ) ) ; 

return (temp) ; 

case 20: 


>Phi dot)*L- 


>Phi dot) *L- 


>Phi dot) *L- 


,->Phi dot) *L- 


,->Phi dot) *T>- 


->Phi dot)*L- 


,->Phi dot) *L- 


->Phi dot) *L- 


(L->Phi_dot) *L- 



zo2 (L->Phi_dot) *L- 


>d[i*3+l] ) ) ; 

return (temp) ; 
case 21: 

L->d_weights t i* 3+0 ] 
L->d weights [i*3+l] 


zero one ( range_pbl <L->Phi) *L->dt i *3+0] ) ; 
zero"one (range _pb2 (L->Phr_dot) *L->d(x 3+1] ) 



temp-min ( zero 
>d [i*3+l] ) ) ; 

return (temp) ; 
case 22: 

L->d_weights [i 
L->d_weights [i 
temp— min ( zero 
>d [i*3+l] ) ) ; 

return (temp) 
case 23: 

L->d_weights [i 
L->d_weights [i 
temp-min ( zero 
>d[i*3+l])); 

return (temp) ; 
case 24: 

L->d_weights [i 
L->d_weights [i 
temp -min ( zero 
>d [i*3+l] ) ) ; 

return (temp) ; 
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»_one (range_pbl (L->Phi) *L->d(i*3+0J ) , zero_one (range_pb2 (L->Phi_dot) 



*3+0] - zero_one 
3+1] - zero_one 


(range_pbl (L->Phi) *L->d[i*3+0] ) ; 

_ (range _pm2 (L->Phi_dot) *L->d [i*3+l] ) • 

>_one (range_pbl (L->Phi) *L->d[i*3+0] ) , zero_one ( range J»m2 (L->Phi_dot) *L^ 


*3+0] — zero^one 
*3+1] - zero_one 
>_one (range pbl (L 


(range_pbl (I*->Phi) & L->d[i*3+0] ) ; 

(range_ps2 <L->Phi_dot) *L->d[i*3+l] ) ; 

->Phi) *L->d [1*3+0] ) , zero_ one (range_ps2 (L->Phi 


dot) 


*3+0] - zero_one 
*3+1] - zero one 


(range_pbl (L->Phi) *L->d[i*3+0] ) ; 

_ (range_zo2 (L->Phi_dot) *L->d[i*3+l] ) ; 

>_one (range_pbl (L->Phi) *L->d[i*3+0] ) , zero_one (range zo2 (L->Phi dot) * t. 


) 


) 


double range_calculate_z_array (i, L) 
int i; 

LEARN_CYCLE * L; /*IN : */ 

( 

switch (i) { 
case 0 : 

return ( range_pm3 (L->w [ 0 ] ) ) ; 
case 1 : 

return ( range_pb3 (L->w [ 1 ] ) ) ; 
case 2: 

return (range_pb3 (L->w[2] ) ) ; 
case 3 : 

return (range_pb3 (L->w[3] ) ) ; 
case 4 : 

return (range_ps 3 (L->w[4] ) ) ; 
case 5: 

return (range_pm3 (L->w[5] ) ) ; 
case 6 : 

return ( range_pb3 ( L->w [ 6 ] ) ) ; 
case 7 : 

return (range_pb3 (L->w [7] ) ) ; 
case 8 : 

return ( range_ps3 (L->w [ 8 J ) ) ; 
case 9: 

return (range_ps 3 (L->w [9] ) ) ; 
case 10: 

return ( range_pm3 (L->w [10] ) ) ; 
case 11: 

return (range_pb3 (L->w[ll] ) ) ; 
case 12: 

return (range_zo3 (L->w[12] ) ) ; 
case 13: 

return (range_nb3 (L->w[13] ) ) ; 
case 14: 

return ( range_nm3 (L->w[14] ) ) ; 
case 15: 

return (range_ns 3 (L->w [15] ) ) ; 
case 16: 

return (range_ns 3 (L->w [16] ) ) ; 
case 17 : 



return (range_nb3 (L->w[17] ) ) 
case 18: 

return ( range_nb3 ( L->w [18])) 
case 19: 

return (range_ nm3 (L->w[19] ) ) 
case 20: 

return (range_ns3 (L->w [20] ) ) 
case 21: 

return (range_nb3 (L->w [21] ) ) 
case 22: 

return (range_nb3 (L->w [22 ] ) ) 
case 23: 

return ( range__nb3 ( L->w [23] ) ) 
case 24: 

return ( range_nm3 ( L->w [24])) 




A2. Source code for elevation control 






♦include <stdio . h> 
♦include <math . h> 
♦include <sys /types .h> 


/* EXTERNAL DATA STRUCTURE DEFINITION 


#include " . . /orb_f uzzy/ learn_cycle . h" 

♦define max(x,y) 

( (x 

>- y) 

? x : y ) 

♦define min(x f y) 

( (x 

< y) 

? x : y ) 

♦define Gamma 

0.9 



♦define Beta 

0.2 



♦define Beta_h 

0.05 



♦define Rho 

1.0 



/* 1 July 92 Change 

Rho h 

from 

0.2 to 0.8 

♦define Rho_h 

0.8 



♦define Rhol 

2.0 



♦define Rho_hl 

0.4 



extern double sgn(); 




extern double exp(); 




extern double rnd(); 




extern double sim time(); 




*/ 


learn_elev (L) 

LEARN_CYCLE * L; /*IN : */ 

{ 

int i,il, j, k; 

double elev_match ( ) , elev_calculate_z_array ( ) ; 
double temp ; 


L->x [ 0 ] - (L->Phi+40) /80 ; 

L— >x [ 1 ] - (L->Phi_dot+10) /20 ? 


/* CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 

/* 1 July 1992 - Set Bias to 0 . 0 */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

L->x [2] - 0.0 ; 

L->failure “ 0; 

/* Set up and evaluate the failure criteria */ 

if ( (fabs (L->Phi) > 40) II (fabs (L->Phi_dot) > 10) ) { 

L->failure - -1.; 
if ( L->learn_f lag ) { 

fprintf (stderr, "learn_elev: Failure at %f . \n", sim_time () ) 
) /* end if */ 

) /* end if */ 

/* ccccccccccccccccccccccccccccccccccccccccccccccccc */ 

/* CC 6 November, 1992 CC */ 

/* CC Turn off learning when there is a failure CC */ 

/* CC until there has been no failure for 4 pass CC */ 

/* CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 

X J ->f 0 urth__pass_failure “ L->third _pass_failure ; 

L->third__pass — failure * L->second _j?ass_f ailure ; 

X J ->second_pass_f ailure ■ L->first _pass_failure ; 
L->first_pass_failure - 0 ; 
if (L->failure “ -1) { 

L->f irst_pass_failure - 1 ; 

} /* end if */ 

/* If No Failure in the last 4 passes then turn learning on */ 


if ( ! L->fourth_pass__f ailure && 

! L->third_ pass_failure && 

! L->second — pass_f ailure && 
! L->f irst_ pass_f ailure ) { 

L->learn_f lag « 1 ; 

} /* end if */ 
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/* If Failure first pass, but not second pass, let the learn 
flag stay on, so that the failure is processed. On the next 
pass (where there is a failure second pass, but not third 
pass) then turn learning off until there are no failures for 
four passes */ 

if ( L->second_pass_f ailure && ! (L->third_pass_failure) ) { 

L->learn__f lag - 0 ; 

) /* end if */ 


/* output: state evaluation */ 

for (i - 0; i < 31; i++) 

( 

L->sum » 0.0; 

for ( j - 0; j < 3; j++) 

{ 

L->sum +- L->a[i*3+j] * L->x_old[j]; 

) 

/*ccccccccccccccccccccccccccccccccccccccccccccccccccc 

CC JUNE 12 , 1992 - Change to "learning speed" CC 

CC L->y[i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; CC 

CC JULY 1 , 1992 - Change to "learning speed" 1.0 CC 

CC L->y [i] - 1.0 / (1.0 + exp (-0.1 * L->sum/3 . 0) ) ; CC 

CC JULY 1 , 1992 - Normalize For Rules Only CC 

CC L->y{i] - 1.0 / (1.0 + exp (-1.0 * L->sum/3 . 0) ) ; CC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC * / 
L->y[i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; 

} 

L->sum - 0.0; 

for(i * 0; i < 31; i++) 

{ 

L->sum +- L->c[i] * L->y[i]; 

) 

L->suml - 0.0; 

for ( j - 0;j < 3; j++) 

{ 

L->suml +- L->b [ j ] * L->x_old [ j ] ; 

) 

L->v - L->sum + L->suml; 

/* output: action */ 
for(i “0; i < 31; i++) 

{ 

il-i; 

L->w [i] - elev__match (il, L) ; 

L->zl[i] - elev_calculate_z_array (il, L) ; 

) 

L->numl - 0.0; 

L->denom - 0.0; 

for (i - 0; i < 31; i++) 

{ 


) 


L->numl +- L->w [i] * L->zl[i] * L->f[i] ; 
L->denom +« L->w [i] *L->f [i] ; 



/* JUNE 9, 1992 - CORRECTION ! ! ! */ 

/* Add test for denom very small compared to suml - no rule firing zone */ 
/* L->push - (1000 . 0*f abs (L->denom) <fabs (L->suml) ) ?0 . 0 :L->suml/L->denom; */ 
if (fabs (L->numl) <-0 . 01 II fabs (L->denom) <-0 . 01) { 

L->push -0.0 ; 

} else { 

L->push - L->numl/L->denom ; 

) /* end if */ 

/* output: action confutations completed */ 
for(i “0; i < 31; i++) 

( 


L->sum « 0.0; 

for (j - 0; j < 3; j++) 

L->sum += L->d[i*3+j) * L->x_old[j]; 

/♦ccccccccccccccccccccccccccccccccccccccccccccccccccc 


CC JUNE 12 , 1992 - Change to "learning speed" CC 
CC L->z[i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; CC 
CC JULY 1 , 1992 - Change to "learning speed" CC 
CC L->z(i] - 1.0 / (1.0 + exp (-0.1 * L->sum/3 . 0) ) ; CC 
CC JULY 1 , 1992 - Normalize For Rules Only CC 
CC L->z [i] - 1.0 / (1.0 + exp (-1.0 * L->sum/3 . 0) ) ; CC 


ccccccccccccccccccccccccccccccccccccccccccccccccccc * / 

L->z [i] » 1.0 / (1.0 + exp (-1 . 0 * L->sum) ) ; 

} 

L->sum2 - 0.0; 

L->sum3 - 0.0; 

for(i “ 0; i < 3; i++) 

L->sum2 +- L->e [i] * L->x_old[i]; 

for (i“0;i < 31; i++) 

L->sum3 +- L->f [i] * L->z [i] ; 

/*ccccccccccccccccccccccccccccccccccccccccccccccccccc 


CC JULY 1 , 1992 - Change Normalize of sum4 CC 
CC L->sum4 - L->sum3 + L->sum2; CC 
CC JULY 1 , 1992 - Normalize of sum3 by # rules CC 
CC L->sum4 - L->sum3 / 31.0 + L->sum2 / 3.0 ; 

CC OCT 2 , 1992 - Do NOT normalize for # rules CC 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/ 

L->sum4 “ L— >sum3 + L->sum2 / 3.0 ; 

/*ccccccccccccccccccccccccccccccccccccccccccccccccccc 

CC JUNE 12 , 1992 - Change to "learning speed" CC 

CC L->p - 1.0 / (1.0 + exp (-1.0 * L->sum4) ) ; CC 

CC JULY 1 , 1992 - Change to "learning speed" 1.0 CC 
CC L— >p ” 1.0 / (1.0 + exp (-0.1 * L->sum4/34 . 0) ) ; CC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/ 

L— >p - 1.0 / (1.0 + exp (-1.0 * L->sum4) ) ; 

/* CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 

/* 15 April 1992 - Use temp variable - not push */ 

/* CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 

/* L->push - ( rnd ( ) <- L->p) ? L->push : -L->push; */ 

/* L->unusualness “ (L->push >0) ? 1.0 - L->p : -L->p; */ 

/* CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 

/* 28 Sept 1992 - Set stocastic action to zero out cmd */ 





/* 

/* 



cccccccccccccccccccccccccccccccccccccccccccccccccccc * / 

L->push - ( rnd() <■» ( (L->p+l . 0) /2 . 0) ) ? L->push : -L->push; */ 


/ * ccccccccccccccccccccccccccccccccccccccccccccccccccccc 
28 Sept 1992 - Change definition of unuaualneas to 
correspond with Hamid' s notes . 

L->push - ( rnd ( ) <- ( (L->p+l . 0) /2 . 0) ) ? L->push : 0.0 ; 
L->unuaualneaa * (L->push >0) ? 1.0 - L->p : -L->p; 

*/ 

if ( rnd ( ) <« ( (L->p+l . 0) /2 . 0) ) { 

L->unusualness - 1.0 - L->p ; 

} else { 

L->push -0.0 ; 

L->unusualness - -L->p ; 

) /* end if */ 


/* using new input values and unmodified weights. */ 

/* Use y_new and v_new so not to destroy y and v. */ 

for (i « 0; i < 31; i++) 

{ 

L->sum - 0.0; 

for(j - 0; j < 3; j++) 

{ 

L->sum +- L->a[i*3+j] * L->x[j]; 

} 

/*ccccccccccccccccccccccccccccccccccccccccccccccccccc 

CC JUNE 12 , 1992 - Change to "learning speed" CC 
CC L->y_new[i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; CC 
CC JULY 1 , 1992 - Change to "learning speed" CC 

CC L->y_new{i] * 1.0 / (1.0 + exp (-0.1 * L->sum/3 . 0) ) ;CC 
CC JULY 1 , 1992 - Do Not Normalize sum CC 

CC L->y_new[i] « 1.0 / (1.0 + exp (-1.0 * L->sum/3 . 0) ) ;CC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC * / 
L->y_new(i] * 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; 

} 

L->sum - 0.0; 

L->suml * 0.0; 

L->sum2 « 0.0; 

for(j *0; j < 3; j++) 

L->suml +* L->b [ j ] * L->x [ j ] ; 

for(i - 0; i < 31; i++) 

L->sum2 +- L->c [i] * L->y_new(i]; 

L->sum “ L->suml + L->sum2; 

L->v_new * L->sum; 

/ *CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC* / 
/* 17 September, 1992 */ 

/* This logic depends upon a "crisp" (two valued) failure */ 
/* It has been replaced by "fuzzy" failures */ 

/* action evaluation */ 
if (L->failure) 

L->r_hat * L->failure - L->v; 
else 

L- > r_hat - L->failure + Gamma * L->v_new - L->v; 

/ *CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC* / 
/* L->r_hat - -1.0 * L->failure + 

/* ( 1.0 - L->f ailure ) * Gamma * L->v_new - L->v ; 

/* modification and update to parameters */ 



sgn (L->c [i] ) ; 
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for (i - 0; i < 31; i++) 

< 

L->f actorl - Beta_h * L->r_hat * L->y[i] * (1.0 
L->c[i] +- Beta * L->r_hat * L->y[i]; 


L->y[i]> 


for(j -0; j < 3; j++) 

{ 

L->a [i*3+ j] +« L->factorl * L->x_old[j]; 

} 


ford " 0; i < 3; i++) 

L->b [i] +- Beta * L->r_hat * L->x_old[i] ; 

for(i “ 0; i < 31; i++) 

( 

/*cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
CC 7 July, 1992 : Weight updates of D's by rule firing strength CC 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc*/ 

L->factor2 - Rho_h * L->r_hat * L->z[i] * (1.0 - L->z[i]) * sgn <L->f [i] ) *L->unus 

ualness ; 

for ( j « 0; j < 3; j++) { 
if ( L->learn_f lag ) { 

L->d [i*3+ j ] +- L->factor2 * L->x_old[j] * L->d_weights [i*3+ j ] ; 

} /* end if */ 

} /* end for */ 


} 


{ 


for(i = 0; i < 31; i++) 


/♦CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CC July 2, 1992 : Include Rule Firing Weight in F update CC 

CC L->f[i] +» Rho * L->r_hat * L->unusualness * L->z[i]; CC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/ 
if ( L->learn__f lag ) { 

L->f[i] +« Rho * L->r__hat * L->unusualness * L->z[i] * L->w[i] 

} /* end if */ 

} 

for(i * 0; i < 3; i++) 


{ 


L->e[i] +- Rho * L->r_hat * L->unusualness * L->x_old[i]; 


L->x_old[0] - L“>x [ 0 ] ; 
L->x_old [1] - L->x[l]; 
L~>x old [2] - L->x [2] ; 


} 




double sgn(x) 
double x; 


if (x < 0 . 0) 

return (-1.0); 
else if (x > 0.0) 
return (1.0); 
else 

return (0.0); 


{ 




/* zero_one function returns 0 for negative numbers 

1 for values > 1 
x for values between 0 and 1 */ 

double zero_one(x) 
double x; 

{ 

if (x < 0) return (0,0); 
else if (x > 1) return (1.0); 
else return (x) ; 


/***★*★*★★*★ Membership Function for Elev_Angle *******************************/ 

double elev_nml (x) 
double x; 

{ 

return (rain ( max<( -x-1) /7 , 0.0 ), 1.0 )); 

} 

double eleven si (x) 
double x; 

{ 

if (x <« -1.0) return (min ( max( (x+8.0)/7 , 0.0 ), 1.0 )); 

else return (min( max( -x/1 , 0.0 ), 1.0 )); 

) 

double elev_zol (x) 
double x; 

if (x <* 0) return (min( max( x+1/1 , 0.0 ), 1.0 )); 

else return (min( max ( (-x+16)/16 , 0.0 ), 1.0 )); 

) 

double elev_psl (x) 
double x; 

if (x <« 16) return (min ( max( x/16 , 0.0 ), 1.0 )); 

else return (min( max( (-x+32)/16 , 0.0 ) f 1.0 )); 

} 

do\ible elev_j>ml (x) 
double x; 

return (min ( max(( x-16) /16 , 0.0 ), 1.0 )); 

) 

Elev_Rate Membership Functions ************************ 1 W 

double elev_nm2 (x) 
double x; 

return (min ( max(( -x-.5)/.5 , 0.0 ), 1.0 )); 

} 

double elev__ns2 (x) 
double x; 

if (x <- -.5) return (min< max(( x+l)/.5 , 0.0 ), 1.0 )); 

else return (min( inax{( -x/.5) , 0.0 ) f 1.0 )); 

} 

double elev_zo2 (x) 
double x; 

if (x <- 0) return (min( max({ x+.5)/.5 f 0.0 ) f 1.0 )); 

else return (min { max(< -x+.5)/.5 , 0.0 ), 1.0 )); 



> 

double elev_ps2 (x) 
double x; 

if (x <- .5) return (min { max(x/.5 , 0.0 ), 1.0 )); 

else return (min( max(( -x+l)/.5 , 0.0 ), 1.0 )); 

} 

double elev_pm2 (x) 
double x; 

return (min ( max(( x-.5)/1.5 , 0.0 ), 1.0 )); 

} 

/***************** Defuzzification Process with Accel Membership Functions ************/ 

double elev_nm3 (x) 
double x; 

( 

return (-0.1 - 0.1*x); 

} 

double elev_ns3 (x) 
double x; 

{ 

return (-0 . l*x) ; 

> 

double elev_zo3 (x) 
double x; 

( 

return (0.0) ; 

) 

double elev_ps3 (x) 
double x; 

{ 

return (0 . l*x) ; 

} 

double elev_pm3 (x) 
double x; 
l 

return (0.1 + 0.1*x); 

) 


double elev_match ( i, L) 
int i; 

LEARN_CYCLE *L; /*IN : */ 

< 

double ten?); 
switch (i) { 


CBS6 0 Z 

L->d weights [i*3+0] - zero_one (elev_nml (L->Phi) *L->d[i*3+0j) ; 

L->d weights [i*3+l] - zero_one (elev_zo2 (L->Phr_dot) *L->dU 3+1] ) ; H n^ *t >d 

temp^nin( zero_one (elev_nml (L->Phi) *L->d(i*3+0] ) , zero_one (elev_zo2 (L->Phr_dot) L->d 

[i*3+l] ) ) ; 

return (temp) ; 

° aS L->d weights [i*3+0] - zero_one (elev_nml (L->Phi) *L->d[i*3+0] ) ; 

L->d weights I i*3+l] - zero_one (elev_ns2 (L->Phi_dot) *L->d[i 3+1]) , , 

temp^ain< zero_one (elev_nml (L->Phi) *L->d(i*3+0] ) , zero_one (elev_ns2 (L->Phr_dot) 

[i*3+i] ) ) ; 



return (temp) ; 

case 2 : — 

L->d_weights [i*3+0] - zero__one (elev_nml (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (elev_nm2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp— min { zero_one (elev_nml <L->Phi) *L->d[i*3+0] ) , zero_one (elev_nm2 (L->Phi_dot ) *L-> 
[i*3+l] ) ) ; - 

return (temp) ; 
case 3: 

L->d_weights [i*3+0] - zero_one (elev_nsl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (elev_ps2 (L->Phi_dot ) *L->d[i*3+l] ) ; ~ 

temp-min ( zero__one (elev_nsl (L->Phi) *L->d[i*3+0] ) , zero_one (elev _ps2 (L->Phi_dot) *L->d 

[i*3+l] ) ) ; 

return (temp) ; ^ 

case 4 : 

L->d_weights [1*3+0] - zero^one (elev_nsl (L->Phi) *L->d [i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (elev_zo2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero one (elev_ns 1 (L->Phi ) *L“>d [ i*3+0 ] ) r zero one (elev zo2 ( L—>Phi_dot ) *L->^ 


[1*3+1])); 

return (temp) ; 
case 5: 

L->d_weights [1*3+0] - zero_one (elevens 1 (L->Phi) *L->d [1*3+0] ) ; — 

L->d_weights [i*3+l] - zero_one (elev_ns2 (L->Phi_dot) *L->d [i*3+l] ) ; 

temp-min { zero_one ( elevens 1 (L->Phi) *L->d [i*3+0] ) , zero_one (elev_ns2 (L->Phi_dot) *L-s 
[1*3+1] ) ) ; _ 

return (temp) ; 
case 6 : 


dot) *L->d 


L->d_weights [1*3+0] — zero^one (elev_nsl (L->Phi) *L->d [1*3+0] ) ; 

L->d_weights [1*3+1] - zero_one (elev_nm2 (L->Phi_dot) *L->d[ 1*3+1 ] ) ; ^ 

temp— min ( zero_one (elev_nsl (L->Phi) *L->d [1*3+0] ) , zero^one (elev_nm2 (L->Phi_dot) *L->d 
[1*3+1])); 

return (t emp) ; 

case 7 : ^ 

L->d_weights [1*3+0] - zero_one (elev_zol (L->Phi) *L->d[i*3+0] ) ; 

L->d weights [i*3+l] — zero_one (elev_pm2 (L->Phi_dot) *L->d [1*3+1] ) ; 

temp— min ( zero__one (elev_zol (L->Phi) *L->d [i*3+0] ) , zero_one (elev _pm2 (L->Phi_dot ) *L-' 
[1*3+1])); — 

return (temp) ; 
case 8 : 

L->d_weights [i*3+0] - zero^one (elev_zol (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [ 1 * 3+1 ] - zero_one (elev_ps2 (L->Phi_dot) *L->d[i*3+l] ) 
temp-min ( zero_one (elev_zol (L->Phi) *L->d [1*3+0] ) , zero_one (elev_j5s2 (L->Phi_ 

[1*3+1])); 

return (temp) ; 
case 9: 

L->d_weights [i*3+0] - zero_one (elev_zol (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [ 1*3+1] - zero__one (elev_zo2 (L->Phi_dot) *L->d[i*3+l] ) 
temp-min ( zero_one <elev_zol <L->Phi) *L->d[i*3+0] ) , zero_one <elev_zo2 (L->Phi_ 

[1*3+1])); 

return (temp) ; 
case 10: 

L->d_weights [1*3+0] - zero^one (elev_zol (L->Phi) *L->d[ 1*3+0 ] ) ; 

L->d~weights [1*3+1] - zero_one (elev_ns2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero_one <elev_zol (L->Phi) *L->d[i*3+0] ) , zero_one (elev_ns2 <L->Phi_dot) *L->^ 

[1*3+1])) ; ^ 

return (temp) ; 
case 11: 

L->d_weights [1*3+0] - zero_one (elev_zol (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (elev_nm2 (L->Phi_dot) *L->d[i*3+l] ) 
temp-min( zero_one (elev_zol (L->Phi) *L->d [1*3+0] ) f 
[1*3+1])); 

return (temp) ; 
case 12 : 

L->d weights [i*3+0] - zero_one (elev _psl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (elev_pm2 (L->Phi_dot) *L->d[i*3+l] ) 


dot) *L-: 


zero one (elev nm2 (L->Phi_dot) *L->d 
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temp -min ( zero_one (elev^psl (L->Phi) *L->d[i*3+0] ) , zero^one (elev_jpm2 (L->Phi_dot) *L->d 
[i*3+l])) ; 

return (temp) ; 


case 13: 

L->d_weights [i*3+0] - zero_one (elev__ psl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (elev_j>s2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp— min ( zero^one (elev_psl (L->Phi) *L->d [i*3+0] ) , zero_one (elev__ps2 (L->Phi_dot) *L->d 
U*3+l])); 

return (temp) ; 
case 14: 

L->d_weights [i*3+0] - zero_one (elev_j>sl (L->Phi) *L->d [i*3+0] ) ; 

L->d_weights [i*3+l ) - zero^one (elev_zo2 (L->Phi_dot ) *L->d [i*3+l ] ) ; 

temp— min ( zero_one (elev__psl (L->Phi) *L->d [i*3+0 ] ) f zero_one (elev_zo2 (L->Phi_dot ) *L->d 
[i*3+l])); 

return (temp) ; 
case 15: 

L->d_weights [i*3+0 ] - zero_one (elev_jDsl (L->Phi) *L->d [i*3+0 ] ) ; 

L->d__weights [i*3+l] - zero_one (elev_ns2 (L->Phi_dot ) *L->d [i*3+l] ) ; 

temp— min ( zero one (elev_psl (L->Phi) *L->d [i*3+0] ) , zero_one (elev_ns2 (L->Phi_dot ) *L->d 

[i*3+l] ) ) ; 

return (temp) ; 
case 16: 

L->d_weights [i*3+0] - zero_one (elev_pml (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (elev_pm2 (L->Phi_dot ) *L->d [i*3+l] ) ; 

temp— min ( zero^one (elev_jpml (L->Phi) *L->d[i*3+0] ) , zero_one (elev_j>m2 (L->Phi_dot) *L->d 
[i*3+l] ) ) ; 

return (temp) ; 
case 17 : 

L->d_weights [i*3+0] - zero_one (elev_pml (L->Phi) *L->d [i*3+0 ] ) ; 

L->d_weights [i*3+l] - zero_one (elev_jps2 (L->Phi__dot) *L->d [i*3+l] ) ; 

temp— min ( zero_one (elev_pml (L->Phi) *L->d [i*3+0] ) , zero_ one (elev_j?s2 (L->Phi_dot) *L->d 
[i*3+l] ) ) ; 

return (temp) ; 
case 18: 

L->d_weights [i*3+0] - zero__one (elev_j3ml (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (elev_zo2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero_one (elev_pml (L->Phi) *L->d[i*3+0] ) , zero_one (elev_zo2 (L->Phi_dot) *L->d 
[i*3+l] ) ) ; 

return (temp) ; 

} 

} 

double elev_calculate_z_array (i, L) 
int i ; 

LEARN_CYCLE * L; /*IN : */ 

{ 

switch (i) { 

case 0: 

return (elev_j?m3 (L->w[0] ) ) ; 
case 1: 

return (elev_pm3 (L->w [1] ) ) ; 
case 2 : 

return (elev_ pm3 (L->w [2] ) ) ; 
case 3: 

return (elev_zo3 (L->w [3] ) ) ; 
case 4 : 

return (elev_ps3 (L->w[4] ) ) ; 
case 5: 

return (elevj>s3 (L->w [5] ) ) ; 
case 6 : 

return (elevjpm3 (L->w[6] ) ) ; 
case 7 : 

return (elev_nm3 (L->w [7] ) ) ; 
case 8 : 


leam_eley,c 


return (elevens 3 (L->w [8 ] ) ) ; 
case 9: 

return (elev_zo3 (L->w [ 9] ) ) ; 
case 10: 

return (elev_ps3 (L->w [10] ) ) ; 
case 11: 

return (elev_j>m3 <L->w[ll] ) ) ; 
case 12: 

return (elev_nm3 (L->w [12] ) ) ; 
case 13: 

return (elevens 3 (L->w [ 13] ) ) ; 
case 14: 

return (elevens 3 (L->w [14] ) ) ; 
case 15: 

return (elev_zo3 (L->w [ 15 ] ) ) ; 
case 16: 

return (elev_nm3 (L->w [16] ) ) ; 
case 17 : 

return (elev_nm3 (L->w [17] ) ) ; 
case 18: 

return (elev ns3 (L->w [18] ) ) ; 



A3. Source code for azimuth control 









♦include <math.h> 

♦ include <sys/types .h> 

/* EXTERNAL DATA STRUCTURE DEFINITION */ 


♦include ” . . /orb_fuzzy/learn cycle. h w 

♦define max(x,y) 

( (x 

>- y) 

? X : y ) 

♦define min(x f y) 

( (x 

< y) 

? x : y ) 

♦define Gamma 

0.9 



♦define Beta 

0.2 



♦define Beta_h 

0.05 



♦define Rho 

1.0 



/* 1 July 92 Change 

Rho h 

from 

0.2 to 0.8 

♦define Rho_h 

0.8 



♦define Rhol 

2.0 



♦define Rho hi 

0.4 



extern double sgn(); 




extern double exp(); 




extern double rnd{); 




extern double sim_time(); 



lea rn_a z im ( L ) 




LEARN_CYCLE * L; 

/* 

IN : 

*/ 


{ 

int i,il, j, k; 

double azim_match () , azim_calculate_z_array () ; 
double temp ; 


/* CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 
/* 11 March 1992 - Alter scaling */ 

/* CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 
/* L->x ( 0 ] - L->Phi/20 . 0; */ 

/* L->x [ 1] = L->Phi_dot / 4 . 0 ; */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

/* 08 April 1992 - Alter scaling */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

/* L->x ( 0 ] - L->Phi/10 . 0; */ 

/* L->x[l] - L->Phi_dot/2 . 0 ; */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

L->x[0] = (L->Phi + 4 0 ) / 8 0 . 0 ; 

L->x[l] >= (L->Phi_dot + 10 ) /20 . 0 ; 


/* cccccccccccccccccccccccccccccccccccccccccc */ 

/* 08 April 1992 - Alter Bias */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

/* L->x [2] - 1.00 ; */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

/* 15 April 1992 - Alter Bias to 0.5 */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

/* 1 July 1992 - Alter Bias to 0.0 */ 

/* cccccccccccccccccccccccccccccccccccccccccc */ 

L->x [ 2 ] - 0.0 ; 


L->£ailure - 0; 

/* ccccccccccccccccccccccccccccccccccccccccccccc */ 

/* CC 12 August, 1992 CC */ 

/* CC Change failure criteria from 0.7 & 0.07 CC */ 

/* CC to be 0.5 and 0.05 CC */ 

/* CC 31 August, 1992 CC */ 

/* CC Change back to 0.7 ( 0.07 CC */ 





/* ccccccccccccccccccccccccccccccccccccccccccccc */ 

/* Set up and evaluate the failure criteria */ 

if ( (fabs (L->Phi) > 2.0) || (f abs <L->Phi_dot) > 0.5) ) { 

L->failure - -1.; 
if ( L->learn_f lag ) { 

fprintf (stderr, "learn_azim: Failure at %f . \n", sim_time {) ) 
} /* end if */ 

) /* end if */ 


/* ccccccccccccccccccccccccccccccccccccccccccccccccc */ 

/* CC 6 November, 1992 CC */ 

/* CC Turn off learning when there is a failure CC */ 

/* CC until there has been no failure for 4 pass CC */ 
/* CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC */ 

L->fourth_pass__failure * L->third_pass_f ailure ; 
L->third_pass_f ailure * L->second_j?ass_f ailure ; 
L->second_j3ass_f ailure * L->f irst_pass_f ailure ; 

L— >f irst_pass failure * 0 ; 

if <L->failure “ -1) { 

L->f irst_pass_f ailure * 1 ; 

) /* end if */ 


/* If No Failure in the last 4 passes then turn learning on */ 
if ( ! L->f ourth_ pass_f ailure && 

! L->third_pass_f ailure && 

!L->second_jpass_f ailure && 

!L->first_pass_f ailure ) { 

L->learn_f lag « 1 ; 

) /* end if */ 


/* If Failure first pass, but not second pass, let the learn 
flag stay on, so that the failure is processed. On the next 
pass {where there is a failure second pass, but not third 
pass) then turn learning off until there are no failures for 
four passes */ 

if ( L->second_joass_f ailure && ! (L->third_pass_f ailure) ) { 

L->learn_f lag * 0 ; 

) /* end if */ 


/* output: state evaluation */ 

for (i =0; i < 16; i++) 

{ 

L->sum - 0.0; 

for ( j - 0; j < 3; j++) 

{ 

L->sum +— L->a[i*3+j] * L->x_old[j]; 

} 

L->y[i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; 

) 

L->sum - 0.0; 

for (i - 0; i < 16; i++) 

( 

L->sum +« L->c [i] * L->y [i] ; 

) 

L->suml - 0.0; 

for ( j - 0;j < 3; j++) 

< 

L->suml +« L->b [ j ] * L->x_old[ j] ; 

) 

L->v - L->sum + L->suml; 



/* output: action */ 
ford - 0; i < 16; i++) 


{ 

il-i; 

L->w[i] - azim_match (il, L) ; 

L->zl(i] - azim_calculate_z_array (il, L) ; 

} 

L->numl ” 0.0; 

L->denom - 0.0; 
for(i “0; i < 16; i++) 

{ 

L->numl +- L->w[i] * L->zl[i] * L->f[i] ; 

L->denom +— L->w[i] *L->f [i] ; 


} 

/* JUNE 9, 1992 - CORRECTION ! ! ! */ 

/* Add test for denom very small compared to suml - no rule firing zone */ 
/* L->push — (1000 . 0*fabs (L->denom) <fabs (L->suml) ) ?0 . 0 :L->suml/L->denom; */ 
if (fabs (L->numl) <“0 . 01 II fabs (L->denom) <-0 . 01) { 

L->push - 0.0 ; 

} else { 

L->push - L->numl/L->denom ; 

) /* end if */ 

/* output: action computations completed */ 
for(i “0; i < 16; i++) 

{ 


L->sum - 0.0; 

for (j - 0; j < 3; j++) 

L->sum +- L->d[i*3+j] * L->x_old[j]; 

L->z [ i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; 

1 

L->sum2 - 0.0; 

L->sum3 = 0.0; 

for(i ■= 0; i < 3; i++) 

L->sum2 +- L->e(i] * L->x_old[i]; 

for (i=0;i < 16; i++) 

L->sum3 +- L->f[i] * L->z [i] ; 

/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


CC JULY 1 , 1992 - Change Normalize of sum4 CC 
CC L->sum4 - L->sum3 + L->sum2; CC 
CC JULY 1 , 1992 - Normalize of sum3 by # rules CC 
CC L->sum4 - L->sum3 / 31.0 + L->sum2 / 3.0 ; 

CC OCT 2 , 1992 - Do NOT normalize for # rules CC 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/ 

L->sum4 » L->sum3 + L->sum2 / 3.0 ; 

/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CC JUNE 12 , 1992 - Change to "learning speed" CC 

CC L— >p - 1.0 / (1.0 + exp (-1.0 * L->sum4) ) ; CC 

CC JULY 1 , 1992 - Change to "learning speed" 1.0 CC 
CC L->p - 1.0 / (1.0 + exp (-0.1 * L->sum4/34 . 0) ) ; CC 
ccccccccccccccccccccccccccccccccccccccccccccccccccc * / 
L->p - 1.0 / (1.0 + exp (-1.0 * L->sum4) ) ; 



i$zpW 


/* cccccccccccccccccccccccccccccccccccccccccccccccccccc */ 

/* 15 April 1992 - Use temp variable - not push */ 

/* cccccccccccccccccccccccccccccccccccccccccccccccccccc */ 

/* L->push «■ <rnd() <■* L->p) ? L->push : -L->push; */ 

/* L->unusualness - (L->push >0) ? 1.0 - L->p : -L->p; */ 


/* cccccccccccccccccccccccccccccccccccccccccccccccccccc */ 

/* 28 Sept 1992 - Set stocastic action to zero out cmd */ 

/* cccccccccccccccccccccccccccccccccccccccccccccccccccc */ 

/* L->push * ( rnd() <- ( (L->p+l . 0 ) /2 . 0 ) ) ? L->push : -L->push; */ 


/ * ccccccccccccccccccccccccccccccccccccccccccccccccccccc 

28 Sept 1992 - Change definition of unusualness to 
correspond with Hamid' s notes . 

L->push - ( rnd() <- ( (L->p+l . 0 ) /2 . 0 ) ) ? L->push : 0.0 ; 
L->unusualness * (L->push >0) ? 1.0 - L->p : -L->p; 

*/ 

if ( rnd ( ) <- ( (L->p+l . 0) /2 . 0) ) { 

L->unusualness -1.0 - L->p ; 

} else { 

L->push -0.0 ; 

L->unusualness - -L->p ; 

} /* end if */ 


/* using new input values and unmodified weights. */ 

/* Use y_new and v_new so not to destroy y and v. */ 

for(i - 0; i < 16; i++) 

{ 

L->sum - 0.0; 

for ( j - 0; j < 3; j++) 

{ 

L->sum +- L->a[i*3+j] * L->x[j]; 

} 

L->y_new[i] - 1.0 / (1.0 + exp (-1.0 * L->sum) ) ; 

) 

L->sum - 0.0; 

L->suml - 0.0; 

L->sum2 « 0.0; 

for(j - 0; j < 3; j++) 

L->suml +- L->b[j] * L->x [ j ] ; 

for(i - 0; i < 16; i++) 

L->sum2 +« L->c [i] * L->y_new[i]; 

L->sum — L->suml + L->sum2 ; 

L->v_new — L->sum; 

/ *CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC^/ 
/* 17 September, 1992 */ 

/* This logic depends upon a "crisp" (two valued) failure */ 
/* It has been replaced by "fuzzy" failures */ 

/* action evaluation */ 
if (L->failure) 

L->r_hat - L->failure - L->v; 
else 

L->r_hat * L->failure + Gamma * L->v_new - L->v; 

/ *CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC* / 
/* L->rJhat - -1.0 * L->failure + 

/★ ( 1.0 - L->failure ) * Gamma * L->v_new - L->v ; 


/* modification and update to parameters */ 
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for(i - 0; i < 16; i++) 

L->factorl - Beta_h * L->r_hat * L->y[i] * <1.0 - L->y[i]) * sgn (L->c [i] ) ; 

L->c[i] +- Beta * L->r_hat * L->y[i]; 

for(j "0; j < 3; j++) 

< 

L->a [i*3+ j] +- L->factorl * L->x_old[j]; 

} 

} 


for(i - 0; i < 3; i++) 

L->b[i] +“ Beta * L->r_hat * L->x_old[i] ; 

for(i - 0; i < 16; i++) 

< 

/*cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
CC 7 July, 1992 : Weight updates of D' s by rule firing strength CC 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc*/ 

L->factor2 - Rho_h * L->r_hat * L->z[i] * (1.0 - L->z[i]) * sgn (L->f [i] ) *L->unus 

ualness ; 

for(j - 0; j < 3; j++) < 
if ( L->learn_f lag ) { 

L->d[i*3+j] +- L->factor2 * L->x_old[j] * L->d_weights [i*3+ j] ; 

} /* end if * / 

} /* end for */ 


} 


{ 


for(i “ 0; i < 16; i++) 


/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CC July 2, 1992 : Include Rule Firing Weight in F update CC 

CC L->f[i] +- Rho * L->r_hat * L->unusualness * L->z[i]; CC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCV 

if ( L->learn_f lag ) < 

L->f[i] +- Rho * L->r_hat * L->unusualness * L->z[i] * L->w[i] 
) /* end if */ 

} 

for (i - 0; i < 3; i++) 


{ 


L->e[i] +“ Rho * L->r_hat * L->unusualness * L->x_old[i] ; 


L->x_old[0] - L->x [0] ; 
L->x_old[l] - L->x[l] ; 
L->x old <2] - L->x [2] ; 


} 






*/ 


double sgn(x) 
double x; 

{ 

if (x < 0.0) 

return (-1.0); 
else if (x > 0.0) 
return (1.0) ; 
else 

return (0.0); 



} 


/* zero_one function returns 0 for negative numbers 

1 for values > 1 
x for values between 0 and 1 */ 

double zero_one ( x) 
double x; 

{ 

if (x < 0) return (0.0); 
else if (x > 1) return (1.0); 
else return (x) ; 


Membership Function for Azim ******************★*★★*★★★**★★★★★/ 

double azim__nml (x) 
double x; 

{ 

return (lain (max{ (-x-2)/2, 0.0) , 1.0)); 

} 

double azim_nsl(x) 
double x; 

{ 

if <x <- -2.0) return (min ( max ( (x+4.0)/2.0 r 0.0 ), 1.0 )); 

else return (min( max( -x/2 , 0.0 ), 1.0 )); 

) 

double azim_zol (x) 
double x ; 

{ 

if (x <« 0) return (min ( max ( (x+2)/2 , 0.0 ), 1.0 )); 

else return (min ( max( (-x+2)/2 , 0.0 ), 1.0 )); 

} 

double azim_psl(x) 
double x; 

{ 

if (x <- 2) return (min ( max( x/2 , 0.0 ), 1.0 )); 

else return (min ( max( (-x+4)/2 f 0.0 ), 1.0 )); 

) 

double a z im_pml (x) 
double x; 

{ 

return (min ( max(( x-2)/2 / 0.0 ) f 1.0 )); 

} 


Azim_Rate Membership Functions 




double azim_nm2 (x) 
double x; 

{ 

return (min ( max(( -x-0.5)/0.5 f 0.0 ) f 1.0 )); 

) 

double azim_ns2 (x) 
double x; 

{ 

if (x <* -.5) return (min( max(( x+l)/.5 , 0.0 ), 1.0 )); 

else return (min ( max ( -x/ .5 , 0.0 ), 1.0 )); 

) 

double azim_zo2 (x) 
double x; 

{ 

if (x <* 0) return (min( max({ x+.5)/.5 f 0.0 ), 1.0 )); 




else return (min { max ( < -x+.5)/.5 , 0.0 ), 1.0 )); 


} 

double azim_ps2 (x) 
double x; 

{ 

if (x <* .5) return (min ( max(x/.5 , 0.0 ), 1.0 )); 

else return (min ( max(( -x+l)/.5 , 0.0 ), 1.0 )); 

} 

double azim_pm2 (x) 
double x; 

< 

return (min ( max(( x-.5)/1.5 , 0.0 ), 1.0 )); 

) 


/******★***★*★*★*★ Defuzzification Process with Accel Membership Functions 


***★***★★**★ j 


double azim_nm3 (x) 
double x; 

{ 

return (-0.1 - 0.1 * x) ; 

} 

double azim_jns3 (x) 
double x; 

{ 

return (-0 . l*x) ; 

} 

double azim_zo3 (x) 
double x; 

{ 

return (0.0) ; 

} 

double azim_ps3 (x) 
double x; 

{ 

return (0 . l*x) ; 

) 

double a z im_pm3 (x) 
double x; 

{ 

return (0.1 + 0.1*x); 

> 


double azim__match (i, L) 
int i; 

LEARN_CYCLE *L; /*IN : */ 

{ 

double temp; 
switch (i) { 

case 0 : 

L->d_weights [i*3+0] - zero_one (azim_nml (L->Phi) *L->d [i*3+0] ) ; 

L->d”weights [i*3+l] * zero_one (azim_zo2 (L->Phi_dot ) *L->d [i*3+l] ) ; 

temp«min( zero^one (azim_nml (L->Phi) *L->d[i*3+0] ) , zero_one (azim_zo2 (L->Phi_dot) *L->d 

[i*3+l] ) ) ; 

return (temp) ; 
case 1 : 

L->d_weights [i*3+0] - zero_one (azim_nml (L->Phi) *L->d[i*3+0] ) ; 

L ->d_weights [i*3+l] - zero_one (azim_ns2 (L->Phi_dot) *L->d [i*3+l] ) ; 

temp -min ( zero_one (azim__nml (L->Phi) *L->d[i*3+0] ) , zero_one (azim_ns2 (L->Phi_dot) *L->d 


,v'; 
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[i*3+l] ) ) ; 


return (temp) ; _ 

case 2 : 

L-*>d__weights [ i*3+0 ] * zero_one (azim__nml (L->Phi) *L->d [ i*3+0 ] ) ; 

L->d_weights [i*3+l] - zero^one (azim_mti2 (L->Phi_dot ) *L->d [i*3+l] ) ; 

temp -min ( zero_one ( a z im_nml (I»->Phi) *L->d [i*3+0 ] ) , zero_one ( a z im_nm2 (L->Phi_dot ) 

[i*3+l] ) ) ; 

return (temp) ; 
case 3: 

L->d_weights [i*3+0] - zero_one (azim_nsl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (azim_zo2 (L->Phi_ dot) *L->d[i*3+l] ) ; 

temp-min ( zero^one (azim_nsl (L->Phi) *L->d [i*3+0] ) , zero_one (azim_zo2 (L->Phi_dot) *L->< 
[i*3+l] ) ) ; _ 

return (temp) ; 
case 4 : 

L->d_weights [i*3+0] - zero_one (azim_nsl (L->Phi) *L->d (i*3+0 ] ) ; 

L->d_weights [i*3+l] - zero^one (azim_ns2 (L->Phi_dot) *L->d [ i*3+l] ) ; _ 

temp-min ( zero_one (azim_nsl (L->Phi) *L->d[i*3+0] ) , zero_one (azim_ns2 (L->Phi_dot) *L->d 
[i*3+l])); 

return (temp) ; 

case 5 : — 

L->d_weights [1*3+0] - zero_one (azim_nsl (L->Phi) *L->d[i*3+0] ) ; 

L->d_weights [i*3+l] - zero_one (azim_nm2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp— min ( zero_one (azim_nsl (L->Phi) *L->d [i*3+0] ) , zero_one (azim_nm2 (L->Phi_dot) *L->< 
[i*3+l] ) ) ; 

return (temp) ; 
case 6 : 

L->d_weights [i*3+0] - zero_one (azim_zol (L->Phi) *L->d[i*3+0] ) ; ^ 

L->d — weights [i*3+l] - zero_one (azim_jpm2 (L->Phi_ dot) *L->d[i*3+l] ) ; 

temp— min ( zero_one (azim__zol (L->Phi) *L->d [i*3+0] ) , zero__one (azim_pm2 (L->Phi_dot) *L->d 
[i*3+l] ) ) ; 

return (temp) ; w 

case 7 : 

L->d__ weights [i*3+0] - zero_one (azim_zol (L->Phi) *L->d [i*3+0 ] ) ; 

L->d_weights [i*3+l] - zero_one (aziinjps2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero_one (azim_zol (L->Phi) *L->d[i*3+0] ) , zero_one (azim _ps2 (L->Phi_dot) *L->w 
[i*3+l] ) ) ; 

return (temp) ; 
case 8 : 

L->d_weights [i*3+0] - zero_one (azim_zol (L->Phi) *L->d[i*3+0] ) ; 

L->d~weights [i*3+l] - zero_one (azim_zo2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero_one (azim__zol (L->Phi) *L->d[i*3+0] ) , zero_one (azim_zo2 (L->Phi_dot) *L->< 
[i*3+l] ) ) ; _ 

return (temp) ; 
case 9 : 

L->d_weights [i*3+0] - zero_one (azim__zol (L->Phi) *L->d [i*3+0 ] ) ; 

L->d_weights [i*3+l] - zero_one (azim_ns2 (L->Phi_dot) *L->d[i*3+l] ) ; _ 

temp-min { zero_one (azim_zol (L->Phi) *L->d[i*3+0] ) , zero_one <azim_ns2 <L->Phi_dot) *L->d 

[i*3+l] ) ) ; 

return (temp) ; 

case 10 : “ 

L->d_weights [i*3+0] - zero^one (azim_zol (L->Phi) *L->d[i*3+0 ] ) ; 

L->d_weights [i*3+l] - zero^one (azim_nm2 (L->Phi_dot ) *L->d [i*3+l] ) ; 

temp-min ( zero_one (azim_zol <L->Phi) *L->d[i*3+0] ) , zero_one (azim_nm2 (L->Phi_dot) *L-> 
[i*3+l] ) ) ; ~ 

return (temp) ; 
case 11: 

L->d_weights [i*3+0] - zero_one (azimjpsl (L->Phi) *L->d[i*3+0] ) ; _ 

L->d weights [i*3+l] - zero_one (aziirrpm2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp-min ( zero_one (azim_j>sl <L->Phi) *L->d[i*3+0] ) , zero_one (azim_pm2 (L->Phi_dot) *L->d 
[i*3+l] ) ) ; 

return (temp) ; — 

case 12: 

L->d_weights [i*3+0] - zero_one (azimjpsl (L->Phi) *L->d[i*3+0] ) ; 
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L->d weights [i*3+l] - zero_one (azim _ps2 (L->Phi_dot ) *L->d [ i*3+l ] ) ; 

temp=min ( zero_one (azim_j5sl (L->Phi) *L->d(i*3+0] ) , zero_one (azim_jps2 (L->Phi_dot) 

[i*3+l]>>; 

return (temp) ; 
case 13: 

L->d_weights [i*3+0] - zero_one (azim _psl (L->Phi) *L->d[i*3+0] ) ; 

L->d""~ weights [ i*3+l ] - zero_ one (azim_zo2 (L->Phi_dot ) *L->d [i*3+l] ) ; 

temp— min ( zero_one <azim_psl (L->Phi) *L->d[i*3+0] ) , zero_one (azim_zo2 (L->Phi_dot) 

[i*3+l] ) ) ; 

return (temp) ; 
case 14: 

L-><i_weights [i*3+0] - zero_one (azim_pml (L->Phi) *L->d(i*3+0] ) ; 

L->d _ weights [i*3+l] - zero_one (azim_pm2 <L->Phi_dot) *L->d [i*3+l] ) ; 

temp=min ( zero_one (azim_pml (L->Phi) *L->d[i*3+0] ) , zero_one (azim_pm2 (L->Phi_dot) 

[i*3+l] ) ) ; 

return (temp) ; 
case 15 : 

L->d_weights [i*3+0] - zero_one (azim_pml (L->Phi) *L->d[r*3+0 ] ) ; 

L->d weights [i*3+l] = zero_one (azim_ps2 (L->Phi_dot) *L->d(i*3+l] ) ; 

temp^min ( zero_one (azim_pml (L->Phi) *L->d[i*3+0] ) , zero_one <azim_ps2 (L->Phi_dot) 

[i*3+l] ) ) ; 

return (temp) ; 
case 16: 

L->d_weights [i*3+0] - zero^one (azimjpml (L->Phi) *L->d [i*3+0 ] ) ; 

L->d weights [i*3+l] - zero_one (azim_zo2 (L->Phi_dot) *L->d[i*3+l] ) ; 

temp— min ( zero_ one (azim_pml (L->Phi) *L->d[i*3+0] ) f zero_ one (azim_zo2 (L->Phi_dot) 

[i*3+l] ) ) ; 

return (temp) ; 

> 

) 



*L->d 


*L->d 


*L->d 


*L->d 


*L->d 


double azim_calculate__z_array (i, L) 
int i ; 

LEARN_ CYCLE * L ; /* IN : */ 

{ 

switch (i) { 

case 0 : 

return (azim_j?m3 (L->w [ 0 ] ) ) ; 
case 1 : 

return (azim_j3m3 (L->w [ 1 ] ) ) ; 
case 2 : 

return (azim_jpm3 (L->w [2 ] ) ) ; 
case 3 : 

return (azim_j?s3 (L->w [3] ) ) ; 
case 4 : 

return (azimj>s3 (L->w [4] ) ) ; 
case 5 : 

return (azim_pm3 (L->w [ 5] ) ) ; 
case 6 : 

return (azim_nm3 (L->w [ 6 ] ) ) ; 
case 7 : 

return (azim_ns3 (L->w [7 ] ) ) ; 
case 8 : 

return (azim_zo3 (L->w [ 8 ] ) ) ; 
case 9 : 

return (azim_j>s3 (L->w [9] ) ) ; 
case 10: 

return (azim_pm3 (L->w [ 10 ] ) ) ; 
case 11: 

return (azim_nm3 (L->w [ 11 ] ) ) ; 
case 12 : 

return (azim_ns3 (L->w [ 12 ] ) ) ; 
case 13: 

return (azim ns3 (L->w [13] ) ) ; 
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case 14: 

return (azim_nm3 (L->w [ 14] ) ) ; 
case 15: 

return (azim_nm3 (L->w [15] ) ) ; 
case 16: 

return (aziin_nni3 (L->w [16])); 



} 


) 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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RUN: V Bar Approach 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


failure vs TIME 

RUN: V Bar Approach 
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SIMULATION APPLICATION: ARIC Hansladonal Coroolta Smokoon 
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SIMULATION APPLICATION: ARIC TVanslational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


f[15] vs TIME 

RUN: V Bar Approach 



TIME (sec) 

MODULE: ORBl TERJm_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


18 


Tbe Nov 17 1992 02:06:07 PM 




0 200 400 600 800 1000 1200 1400 1600 1800 

TIME (sec) 

MODULE: ORBITERim_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 19 Tic Nov 17 1992 02:06:07 PM 












TIME (sec) 

MODULE: ORB ITER-lm_el ev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


22 


Tbe Nov 17 1992 02:06:07 PM 


SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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RUN: V Bar Approach 



TIME (sec) 

MODULE: ORBlTERim_elev 

DATA SAMPLING FREQUENCY: 0.200 Hz 


ORBITAL OPERATIONS SIMULATOR 


23 


Tie Nov 17 1992 02:06:07 PM 


SIMULATION APPLICATION: ARIC Translational ControUcr Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational 
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SIMULATION APPLICATION: ARIC TVanslational Controller Simulation 


d[27] vs TIME 

RUN: V Bar Approach 
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SIMULATION APPLICATION: ARIC TVanslational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


push vs TIME 

RUN: V Bar Approach 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


x[l] vs TIME 

RUN: V Bar Approach 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 


d[39] vs TIME 

RUN: V Bar Approach 
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SIMULATION APPLICATION: ARIC Translational Controller Simulation 
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